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I. l7itroduction. 



In the year 1839 Gauss published his celebrated Memoir on Terrestrial Magnetism, 
in which the potential on the Earth's surface was calculated to 26 terms of a series of 
surface harmonics. It was proved in this Memoir that, if the horizontal components 
of magnetic force were known all over the Earth, the surface potential could be derived 
without the help of the vertical forces, and it is well known now how these latter can 
be used to separate the terms of the potential which depend on internal from, those 
which depend on external sources. Nevertheless Gauss made use of the vertical 
forces in his calculations of the surface potential in order to ensure a greater degree 
of accuracy. He assumed for this purpose that magnetic matter was distributed 
through the interior of the Earth, and mentions the fair agreement between calculated 
and observed facts as a justification of his assumption. In the latter part of the 
Memoir it was suggested that the same method should be employed in the investiga- 
tion of the regular and secular variations. 

The use of harmonic analysis to separate internal from external causes has never 
been put to a practical test, but it seems to me to be especially well adapted to 
enquiries on the causes of the periodic oscillations of the magnetic needle. 

If the magnetic effects can be fairly represented by a single term in the series of 
harmonics as far as the horizontal forces are concerned, there should be no doubt as to 
the location of the disturbing cause, for the vertical force should be in the opposite 
direction if the origin is outside from what it should be if the origin is inside the 
Earth. As the expression for the potential contains in one case the distance from the 
Earth s centre in the numerator, in the other case in the denominator, and as the vertical 
force depends on the differential coefiicient with regard to distance from the Earth s 
centre, each single term in the series is of opposite sign according to the location of the 
cause ; but what is true for each single term need not be true for the sum of the series. 
By a curious combination of terms the vertical forces might possibly be of the same 
sign, on whichever of the two hypotheses it is calculated. In any case, however, the 
differences between the two results will be of the same order of magnitude as the 
vertical force itself. If it is then a question simply of deciding whether the cause is 
outside or inside, without taking into account a possible combination of both causes, 
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the result should not be doubtful, even if we have only an approximate knowledge of 
the vertical forces. 

Two years ago I showed that the leading features of the horizontal components for 
diurnal variation could be approximately represented by the surface harmonic of the 
second degree and first type^ and that the vertical variation agreed in direction and 
phase with the calculation on the assumption that the seat of the force is outside the 
Earth. The agreement seemed to me to be sufficiently good to justify the conclusion 
that the greater part of the variation is due to causes outside the Earth's surface. 
Nevertheless, it seemed advisable to enter more fully into the matter, as in the first 
approximate treatment of the subject a number of important questions had to be left 
untouched. I now publish the results of an investigation which has been carried as 
far as the observations at my disposal have allowed me to do. My original conclusions 
have been fully confirmed, and some further information has been obtained which I 
believe to be of importance. The results of the calculation point not only to an 
external source, but to an additional internal source, standing in fixed relationship to 
the external cause. This we might have expected. A varying potential due to 
external causes must be accompanied by currents induced in the Earth's body, which^ 
in turn, must affect the magnetic n€jedle. The phase of these currents and their 
magnitude lead us to form definite conclusions on the average conducting power of 
the Earth, and it will be seen that there is strong evidence that the average con- 
ductivity is ver}'- small near the surface, but must be greater further down. In this 
part of the investigation I had much assistance, from my colleague. Professor Lamb. 

I hope that the results obtained in this paper may induce the heads of magnetie 
observatories to consider the suggestions which I have made at the end of it, as their 
adoption would very materially assist further investigations, 

I had, in the first place, to fix on a year for which we possess as complete magnetic 
records as possible. The phase of the variation of horizontal force changes sign in a 
latitude not far removed from that of Lisbon, and it seemed to me, therefore^ essential 
that the excellent observations there made by Sen. Jos. Capello should be made use 
of The observations are published as far as 1872, and I had to take a year therefore 
anterior to this. It seemed al^o desirable to make use of the St. Petersburg observa- 
tions, as it is the most northerly station for which we have records extending over a 
period of years, and as Mr. H, Wild'b well-known skill gives special value to the 
observations made under his direction. The obseravtions (continued since 1878 in 
Paulowsk) were interrupted in 1871 and 1872, and we have to go back, therefore, as 
far as 1870 if we want to utilise the St. Petersburg and Lisbon observations simul- 
taneously. As far as the horizontal components are concerned, we also possess good 
records of 1870 at Greenwich and Bombay. Four stations are sufficient to find to 
the necessary accuracy the potential on the surface of the Earth, but it would be of 
advantage if in future similar investigations a greater number of stations could he 
utilised. 
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The observations are published in very different form by the various observatories. 
The units of force at Bombay and St. Petersburg are the Gaussian unit millimeter- 
milligram-second. At Lisbon it is the foot-grain-second. At Greenwich the variations 
are given in terms of the whole vertical and horizontal force. At Bombay, moreover, 
the observations are given not for a certain hour, but for a time which varies with 
different instruments between 12 and 1 9 minutes past each hour, For this there was 
some reason originally, but ^t present it would be far better if Bombay w^ould adopt 
the practiQe of other observatories. The daily variations had, in the first place, to be all 
reduced to C.G.S. units, and, further, instead of variations in declination and horizontal 
force, we had to find the components of the periodic force towards the geographical 
North and West, I am much indebted to Mr. Wm. Ellis, of the Greenwich Observatory, 
for the help he has given me in the reduction of the observations to a form which I 
could use in my calculations. ^ good many of the computations were done under his 
direct superintendence, and much time and trouble was saved me in consequence. 

The daily variation of declination and horizontal force was expressed in the form 

ary cost + hi aiii t + c?^ cos 2t + h^ sin 2t + a^ cos 3t -{- &3 sin M + a^ cos U -f b^ cos 4^, 

where t represents astronomical time. The summer months, April to September, were 
treated separately from the winter nionths, October to M^roh. The unit of force, for 
convenience's saj^e, w^s taken as 10"^ C.G.S. 

Tables I, and II. give the cpefScients which were ealeulated according to a well- 
known method from the original observations. 

Tables III, and IV. give the same co^flScients reduced to forcps directed to the 
geographical North §^nd West, instead of to the magnetic North and West. 

Before showing how, with the help of these coeflacients, the sqrfacp potential can be 
calculated, I mui^t deduce a few formulae which v^ill be used hereafter. 

Tabids I. — Force to Magnetic West. 
Coefficients in the expansion 

a^ cos t + ^1 sin t+a.2 ^^^ ^^ + \ ^i^ ^^ + % ^^^ 3^ + \ sin 3^ + a^ cos 4zt + 6^ sin At. 
The unit of for!0o is 1 C.G.S. X 10"^^; t being astrononqiical time. 
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Table II.— Force to Magnetic North. 

Coefficients in the expansion 

a I cos I + &i sin i + % cos 2t + 5^ sin 2t 4- % cos 3^ + h^ sin 3^ + % cos 4^ + h^ sin 4^. 

The unit offeree is 1 C.G.S. X 10~®; t being astronomical time. 
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Table III.— Force to Geographical West 

Coefficients in the expansion 

% cos ]{ + &i sin « + (I3 cos M + \ sin 2t + % cos 3^ + \ sin M + % cos U + h^ sin U. 

The unit offeree is 1 C.G.S. X 10*"^ ; ^ being astronomical time. 
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Table IV. — Force to Geographical North. 

Coefficients in the expansion 
Oi cos t +h^Bmt -^ %cos 2t + h sin 2^ + % cos St + h^ sin 3t + a^_ cos it + h^ ein it 
Tlie unit offeree is 1 O.G.S. X 10""^ ; t being astronomical time^ 
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n. Some Fonnulce useful in the Analysis by Spherical Harmonics, 

In the expansion of mathematical expressions into spherical harmonics, it will often 
occur that we have to express powers of the cosine of an angle, or a cosine of some 
multiple of an angle, in terms of the differential coefficients of zonal harmonics with 
respect to the cosine of the argument. The necessary equations for the powers of 
cosines can be easily obtained by differentiation of the well-known fomiulse, giving 
the powers of a cosine in terms of zonal harmonics. But I think it will be useful here 
to give the general equations which I have had to use in expressing Gosmtt in terms 
of d^Vi/dfiP, where p is any given number, /x = cos u and P,- the zonal harmonic of 
degree i. 

We may start from the expression 



cos mu 



where 



■^m-^.m I ^m — 2'^^i — 9> J • * • -^i^i i * 



(1). 



^ "*" ^ {m - (i + l)}{m -{i-~l)} ... (m ~T) (m +T) . . . {m +Ji + 1)}' 



li' ni and i be even, and 






('9 ■ -J_ -i\ ^''^ — ( i — 2)} {jft — (f — 4)1 ■■ . (m — 1) (m + 1) ■ . ■ {m + {i — 2 )} 
^ '^'^ ' (m -if+ l)}{m - {i -1)} . . . (m -~2) (m + 2) . . . {m + (i -t 1)}' 



if m and i be odd. 
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Both expressions are included in the general one 

^ ^ ^ {ni -^(i + l)}{m ^(i^l)},,, {,n + (^ + 1)} 

where successive factors of both numerator and denominator increase by 2.'"^ 
If in formula (l) we substitute for any term P| 

//p. fJT>. 

Hl±l _ "tllizl ^ /2i JL I) p. (o\ 

dfi dfi i^^e. -r M^^ ••••••. . {^), 

we obtain an expression for coBm0 in terms of the first differential coefficients of the 
zonal harmonics ; and^ if in the formula so obtained we successively apply the trans- 
formation 

WB finally obtain the required expression in tdrms of any required differential 
coefficient. 

I have obtained in tbis %ay the equation 

where d^Pi is written shortly for d^Fi/d^jf, and 

A,=^(2i + ^m^'^^ (4) 

^ ^ ^ {m-(^ + p + l)Hm- (^+j?-l)}.. .{m + (i+p + 1)} * ^^' 

except for the last ^errb. of the series, which will be^given presently 1(5). If p = 0, 

this expression agrees with th^ one previously given, and, as I shall proceed to show, 

if it is true for any valu^ p, it will also be^rue forthe valub_p +^1. Assume, then, 

tb e eq nation (3) to hold . 

The delation 

^dP^T^^^ -^dP^'F^^r, = ^i + 1) #P, 

shows thiat the Mtor multiplying # + T^^:i in the expression for ^ 60s m^ will depend 
only on-A^ aild on A^^ ,.?^j in (3). If B^i^i^ is this fector, we lia;ve 



B 



-*rX.5 JLJL'7 



% -\- Jj 



^' + ^> 2i+'l 2i + 5 



5 



'* In the very useful book on SSpherical Harmonics/ bj Fieeies, the factor m does not occur in the 
general expression for Ai (page 33) ; but, from the deduction of the formula, it is clear that the factor m 
must be taketi twice wheu it occurs in the numerator, and not at all when it occurs in the denominator* 
In using the equation I was at first led into error hj the ambiguity, and hence I believe the expression 
given above to be clearer. 
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or, by substituting A^ from (4) and A^-^^ from tbe corresponding equation, 



B,' . 1 = m 



(m — (i — p — 2)} {m —■ (i — p — 4:)} ... {m + (^ — ^ — 2) } 



H-{-l 



{m — (^ + J? + 1) } {w^ — {i — p 



The square bracket reduces to 

{i — pY — (i -}• p + sy 

{m — (^ + p + 3) } {m + (i + p + o)} 

so that 

cos m0 



- 1) } . . . {m 4- (i + J? + 1) } 

{m — (i -{■ p + S)} {m + (i + p + S)}_ 



(2i + 3) (2p + 3) 



(m — (i 4- J? + 3) } {m + (i + p + S)} 



X . o . , , . 

where 



(2p + 3) 



A r/^ + ip 4- A. rJp + ^V- 



A^+i = m(2i+ 3) 



{m —• (i — p — 2)} {m — (i — p — 4:)} ... {m + (i — p — 2) } 
{m — (i + j^ + 3) } {m — (-i + i^ + 1) } . . . {m + (i + ^ + 3) } 



which agrees with (4) if p + 1 is written for p, and ^ + 1 for i. 

The end terms require special consideration. 

As only the even or only the odd zonal harmonics enter into any one series, the 
difference between m + P ^^d i in the equation (3) must be even ; hence, the 
numerators in the fractional expression of (4) consist always of even numbers, what- 
ever values m, i, or p may have. 

If (m + p) he odd^ the zonal harmonics will all be odd, and the last term of the 
series will depend on d^Vp + i if p be even, and d^P^ if p be odd, for the differential 
coefficients of the zonal harmonics of lower degree will vanish. The expression (4) in 
neither case gives the correct factor, and we must substitute for it in both cases 



(2^ + 1) 



(m — (i -f- p + 1)} {m — (i -j- p — 1)} ... {m + (i -}- p -{■ 1)} 



(5). 



If in (4) i is put equal to (p + 2), the first and last term will be equal to m; in that 
case the m is only taken once. 

The expression (5) is easily proved, and shown to hold also if (m + p) be even. 

The first term of the series is included in the general expression (4). 

I have deduced with the help of the equation (3) the following relatiojjs, which will 
be used in this paper : — 



cos u 






(A), 



cos 2u 



JVll-'OOOljJv-A. Jki^,"""""^ • 



4 ^ 
X 5 



dfi 



Ts";^ . . . . 



• ■ (B), 



3 p 
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COS 3u 



_ 8.^ 4 3 

3 5 7.. "" 7 



dfl dfl 



• * 



(C), 



JT> ^jjy ^p 

/7P /7P /7P 

C0SDt^^e93^^ 77 ^^_^ i- e3 ^^ . , » {S^), 

^^ca A'?/ — l-lt_. ^?I 1 2 8 ^^ } „4„ ^?B I 1 ^Jj /T^^\ . 

COS ()^^ ^ 300-3 ^^ - 455 ,^^ + 55 ^^ + 1 5 ^^ l^^ j^ 

COB tl — -Yo^ ^ ..,......,....,, \^)^ 

^3p ^73 p 

p^^ 9., _ _4__ ^^i 7_ ^_t3 /TT^ 

^3p ^73T> 

COS 6U — -sTKI ^/^B "^ 1>"5 7^8 '<»'•••••••'• v-*-;? 






COS 4^^ = ■4 5V4T"^78' "™ TY¥5 "^TTs" « ¥9T 3^1" ' • ' ' • • • •' ' V-^/? 



^Sp ^3p ^/Sp 

COS Dt^-™- 135135^^3 14~85 ^^3 + 429 ^^^8 • ^ - « • • • 1^, 

^^., />,, 512 ^'^9 _i2 8._ ^^^7 I _4__ ^^^, ...^ 7 __ ^Ia / MN 

COS bW - 765765 ^^3 ^ 36465 ^^^b" ^T SB'S ^^3 12 8 7 ^^S * ' ' i^^-^/' 

There is another formula useful in similar investigations, which may find a place 
here, although I have not used it in the final reductions. It is the expression of a 
zonal harmonic in terms of the n^^ differential coefficient of other zonal harmonics. 



(2i + 2n -f 1) (2i + 2?^ - 1) . . . (2^ + 1) ^^ "''^'^ 
n (2i + 2n — 3) 



l{2i-h2n-l){2i + 2n-$), ,.{2t-iy ^ ''- 

n . n-1 (2 i + 2p.— 7) 

1.2" * (27T2^r^^^^^^^ - 3) 



«. 






=*= (2i 4- 1) (2i -"i)T . . (2i zj^;;i;ziij ^"' ^"^' 



the general p*^' term being 

The proof is conducted exactly in the same way as for equation (3). If i is equal to 
or greater than 2n, the series has its full number of terms, viz., ^ + 1, otherwise the 
series breaks off owing to the differentiaJ coefficients vanishing. 
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III. Expansion of Potential in terms of a Series of Surface Harmonics, 

We know from observation that, excepting the Arctic regions, the daily variation of 
the West force is nearly the same along the same circle of latitude. It is this fact 
which renders the present investigation possible, as comparatively few places of 
observation will be necessary to give us a very fair idea of the nature of the oscil- 
lation over a considerable area of the Earth's surface. If at any place we find that the 
daily oscillation of any one element can be expressed in the form 

a^ cos ^ + 6x sin t + a^ cos 2t + &2 ^^^^ 2^ + • - ^ 

when t is recl^oned by local time, we may get the variation at any point of the same 
latitude circle by writing t -\-\ for t, where X is the longitude towards the East from 
some standard meridian and t now is the time of the standard meridian. At the time 
^ = we have then the variation of the force to geographical West in different longi- 
tudes expressed by 

Y =: a^ cos X -f ^1 sin X -\- a<^ cos 2X -f h^ sin 2X. 

The coefficients will be functions of the latitude, and by expressing these functions 
in a series of proper form we may at once obtain an expression for the potential. 

If X is the force to geographical North, Y the force to geographical West, and Z the 
vertical force, reckoned positive upwards, we have, putting u for the colatitude, and X 
for the longitude towards the East, 

^^ dN . dN ^ dV 

A. = — TTj I Sm U ■=■ — r^ ? Z :^ — -— -? 

a WW a dX dr 

a being the Earth^s radius. 

If we can expand Y sin to in terms of surface harmonics our task is accomplished, 
and for this purpose we need only express % sin u, h^ sin u, &g., in a series of tesseral 
harmonics. 

The expansion of a function of an angle in terms of the trigonometrical functions of 
its multiples is so easily carried out by the method of least squares, if the function is 
given for a regular series of submultiples of 27r, that it seemed to me to be the easiest 
method of proceeding to obtain first by interpolation from the observed coefficients of 
Y its values for equidistant circles of latitude. 

To obtain a curve for each of the values a^, b^, a.2, h^, as depending on the 
latitude, we have the following data. The values are directly observed for four 
points in the Northern hemisphere. Taking the potential (as it is observed to be) 
symmetrical on the Northern and Southern hemispheres, and fixing the period to 
which we apply the calculation to be the one for which the mean value of the 
observed summer variations in the Northern hemisphere holds good, we must put, for 
the West force at the same period, in the Southern half of the globe, the observed 
ivinter values of the corresponding Northern latitudes, reversing, however, the sign to 

3 P 2 
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make our expression agree with observation. We have then eight values for the 
different terms of Y through which we might at once proceed to draw a curve ; but, 
making use of the observed variations of Northern force, we can also calculate the 
direction of the tangent of the required curve for the same eight points. This has 
been done as follows. Let Y be expressed as 

t{a^ cos n\ + b^ sin nX), 



and X as 



S(a^i COS n\ + ^n sin nX). 



The existence of a potential implies the relation 

dY sin to dX 



dib 



d\ 



or 



dY dX 



dib 



dX 



cosec u — Y cot % 



from which, by substitution, 
S ( -7- COS nX + ~r- sin nX 



dvj 



du 

= %[{n^n cosec u 



an cot u) GosnX — {na,, cosec u + K cot u) sin nX] ; 



and, therefore, 



du 



db 



n 



chh 



= nj3n cosec u -— a^^ cot u^ 
= — na^ cosec ii + b^ cot u. 



These equations give the rate of change per radian ; to get the rate of change per 
degree of colatitude we have to multiply the differential coefEcients with the circular 
measure of one degree. The following Table gives the rate of change of the West 
force per degree of colatitude, calculated as explained from the North force. 

Table V. 
The unit offeree in this Table is 10"^^ C.G.S., and the unit of latitude is 1 degree. 



dajdu 
dhjjdu 

da^jdu 
dhjdu 
dag/dit 
dhjdu 
dajdu 
dbjdu 



Bombay. 



Summer. 



- 0-201 

- 7-853 

- 1-184 

- 6-079 

- 2-921 

- 2-420 

- 1-125 
+ 0*407 



Winter. 



+ 
-6 
+ 
-4 

- 1 
-2 

- 1 
---1 



714 

326 
258 
754 
968 
277 
768 
117 



Lisbon. 



Summer. 



- 1 210 

+ 0-432 

- 2-087 



+ 



0-948 

1-218 
0-527 
1-084 
0-471 



Winter. 



- 2-561 
+ 0-931 

~ 2-207 
H- 1-684 

- 0-658 
+ 0-580 
+ 0-503 
+ 0-018 



Greenwich. 



Summer. 



+ 0-356 
+ 2-118 
+ 0-150 
+ 5-445 

- 0-416 
4- 0-988 
+ 0-501 

- 1-036 



Winter. 



- 2-793 
+ 2-053 

- 0-330 
4- 3-318 
+ 0-600 
+ 1-564 
+ 0-622 
H- 0-109 



St, Petersburg. 



Summer. 


Winter. ; 


+ 1-951 


3-086 


H- 2-342 


+ 1-107 


0-079 


0-605 


+ 7-185 


+ 2-689 


+ 2-805 


+ 1-650 


0-030 


+ 1-112 


+ 0-329 


+ 640 


f 0-936 


+ 0437 
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Curves were now carefully drawn for each of the eight coefficients, naaking them fit 
in as well as possible with the ordinates and the direction of their tangents, as given 
in Tables III. and V. 

The values of Y were then read off for each 7° '5 of colatitude, and a fresh table 
was formed (Table VI.). 

From this point onwards we have to carry on the calculations separately for each 
type of the variation. 

Table VI. 

Coefiicients in the series Y = aj^ cos t -\- b^ sin t + <^'2 ^^^ 2^ + ^s sin 2t + &c. for 
different degrees of colatitude, the unit of force being C.G.S. X 10""^. 



Colatitude. 


a^. 





a_;. 


h,. 


a,}. 


h- 


( 0) 


K 


o 




+ 10 











(+48) 


( 0) 


7-5 


^ 32 


+ 97 


+ 42 


-f 9 


(+12) 


(+51) 


(+ 1) 


(+ 3) 


15-0 


+ 52 


+ 158 


+ 76 


+ 22 


(+ 24) 


(+53) 


(+ 1) 


(+ 6) 


22-5 


+ 71 


+ 196 


+ 91 


+ 48 


+ 39 


+ 54 


+ 1 


+ 7 


30-0 


+ 85 


+ 215 


+ 90 


+ 101 


+ 56 


+ 54 


+ 2 


+ 11 


37-5 


+ 96 


+ 231 


+ 89 


+ 143 


+ 73 


+ 55 


+ 11 


+ 18 


45-0 


+ 92 


+ 243 


-f 96 


4- 160 


+ 99 


+ 61 


+ 16 


- 3 


62-5 


+ 83 


+ 247 


+ 90 


+ 168 


+ 98 


+ 58 


+ 22 


-21 


60-0 


+ 74 


+ 205 


+ 83 


+ 165 


+ 91 


+ 54 


+ 20 


-24 


67-5 


4- 65 


+ 144 


+ 131 


-f 83 


+ 89 


+ 28 


+ 11 


-20 


75-0 


+ 63 


+ 87 


+ 127 


+ 37 


+ 69 


+ 10 


+ 2 


- 15 


82-5 


+ 58 


+ 58 


+ 117 


+ 15 


+ 46 


- 3 


6 


- 11 


90-0 


+ 37 


+ 42 


-f 82 


+ 7 


+ 25 


- 2 


- 14 


+ 6 


97-5 


+ 10 


+ 21 


+ 15 


+ 1 


+ 3 


4- 1 


20 


+ 18 


105-0 


- 3 


- 7 


- 4 


17 


- 16 


6 


27 


+ 13 


112-5 


+ 3 


- 56 


- 3 


48 


36 


-23 


-40 





120-0 





105 


4 


102 


-47 


45 


-43 


-23 


I27'5 


50 


120 


15 


118 


-53 


-50 


-35 


-26 


135-0 


78 


113 


16 


- 106 


- 54 


-46 


-31 


-25 


142-5 


- 97 


- 100 


+ 10 


- 85 


- 47 


37 


-24 


-24 


150-0 


-121 


- 85 


+ 6 


- 65 


- 38 


30 


-12 


-14 


157-5 ' 


- 140 


75 


+ 2 


- 46 


- 22 


-21 


- 5 


-13 


165-0 ; 


-154 


- 61 


+ 1 


30 


(- 9) 


(-14) 


(- 3) 


(- 9) 


172-5 


- 162 


- 36 





- 14 


(- 2) 


(- 6) 


(- 2) 


(~ 4) 


180-0 


- 165 











( 0) 


( 0) 


( 0) 


( 0) 



The quantities a^, ^^ were expressed in the usual way in terms of the multiples of 
the cosines of the colatitude, and two equations obtained which can be shown to 
represent wdth sufiicient accuracy the force of the first type towards the geographical 
West. 

aj^ z= — . 5 + 1067 cos u — • 50*6 cos 2u =-17 cos Sic — 14'6 cos ^u 

— 8 '9 cos 5u — 2*0 cos 6u, 

b^ = 49*4 + 154'3 cos u — 0*4 cos 2u — 76*0 cos Su — • 18'1 cos 4:U 

•— 31'3 cos 5u — 9'2 cos 6u. 



If, now, cos II, cos 2u, cos 3u, &c., be expressed in terms of dPJdii, dPc^/dfii, dPJdfiy 
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&c., we have an equation for that part of Y which depends on cos X and sin X. 
After multiphcation with sin u^ one side of the equation contains the colatitude in 
form of tesseral harmonics only, and hence we obtain at once the required 
expansion of V^^l The equations are 



Y^^> =:= COS X 



+ sin X 



24-09 



^ + 35-59 ^ ~- 8-45 ^^ + 2-39 ^ 
d^ clfju a/jU ci/jb 



2-41 ~r^ 
cljj, 






48-00 '7^J- + 81-52 '^ + 5-66 "^ 
a/jb a/M dfi 



r/V d P 

7-62 ~P- 1-09^ 
djju dfjb 



578^ 
d/u, 



With the help of 



djx, dfM J 



dY 
dX 



Ya sin u, 



we have, finally 



V^w/a = cos X[48-00Tii + 81 -5217 + 5-661V ~ 7-621V 

- rogT^i - 5-78Tei - 1-571V] 
+ sinX[ - 24-09Tii + 35-59T3I + 8-45131 - 2-39141 

+ 2-411V + 1-641V + 0-3417!] . 



[A]. 



To obtain^ similarly, an expression for V^^^, the series of coejSicients a^, h^ were 
expressed in terms of sin u^ sin 2u, Ac, giving the equations 

a.2 =^ 67 '0 sin u + 61 '3 sin 2it + 7*4 sin Su — 11 '6 sin Au + 15*9 sin 5t^ + 18'6 sin Gu, 
h.^ = 21*5 sin u + 113'2 sin 2u + 77 sin 3ic — 12'6 sin 4.U — 10*8 sin 5u — - 28*2 sin 6?/. 



From the known expansion of cos pu^ in terms of the first differential coefficient of 
the zonal harmonics, we can obtain by differentiation an equation which will at once 
give us V^^^ in the required form. 

Thus, for iDstance, w^e have by Equation D, page 474, 



cos 4it 



315 



d/M 



16. 

4 5 



r/P 



-f ' 3' 



and differentiating with respect to /x, 



4 sin 4,u I sin i^ == 3^ 



#P 



16 _„3 1 3 '^^ ^ 1 



16- ^..3 45 ^^3 -r 35 ^^2 



diM' 



By substituting this and other similar expressions^ we find 
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Y^^'-'sin ti = cos 2X 



64-59 + 9-42 ""ff -~ 0*42 / . ' ~ l'^^ ^ 

dfji'^ ajj.'' dyu 



3 






+ sin 2X 



(PV d^V d^V 

17-6.9+ 15-87 TT + 1-27^^+ 0-68^ 

a/x-^ dyC^ djji'' 



d- P ^2P 

0-40 V„^- 0-80 ^/f 



and from this, by multiplication with sin^ u and integration with respect to X, 

- V^sya = cos 2X [2*95 T/ + 7*94 T^ + 0-63 T/ + 0-34 T/ - 0-20 T/ - 0-41 T/] 

^sin2X[10-76T/ + 4-72T3^~0-2lT/~-073T,H0-29T/ + 0-26T/J. [B]. 

To find V^^^ from the coefficients ag, h.^, ci^, 6^,, the values of these coefficients were all 
divided by sin^ u, for a reason which will appear. It was then found, however, that to 
represent the numbers so obtained satisfactorily, a greater number of terms was 
necessary in the expansion than the observations seemed to justify. I Avas thus led to 
give up calculating as far as these types are concerned the terms on which the differ- 
ence between summer and winter depends, so that, instead of using the coefficients from 
Table VI., the mean of corresponding values on both sides of the equator were taken, 
changing, of course, the sign of the coefficient given for the Southern hemisphere. In 
the next place, it was found that, owing to the smallness of sin u at high latitude, the 
three numbers corresponding to the colatitudes 0"", 7''-5, IS'' were very large, while their 
weight is very small, as they have only been obtained by graphical extrapolation ; 
I have, therefore, discarded these numbers altogether, putting them into brackets in 
Table VI. Bessel has shown how the method of least squares may be applied to 
obtain a trigonometrical series for a succession of values some of which are missino-. 

The equations thus calculated are as follows : — 

ag/sin^ u = 74 cos u -—118 cos Su — • 101 cos 5ii, 
63/sin^ It = 23 cos i^ — ill cos 3t^ — 127 cos 5u, 
ajsin^ u:= 39 sin 2u — 6 sin Au — 9 sin Gu, 
64/sin^ t/^ =1= 31 sin 2u + 38 sin Ati + 4 sin Cru, 

The agreement with the numbers from which the series are obtained is not as good 
as would be desirable, especially for a^ and 63, but the types higher than the second 
Avill probably depend much on local circumstances, and the result would not, in 
my opinion, repay the trouble of taking account of further terms in the above series ; 
cos u, cos 3u, and cos 5u have already been given in terms of d^l\/dii? (Equations 
G, I, L, page 474), and, by differentiation with respect to u, we can also get sin 2ii, 



480 



PROFESSOR A. SCHUSTER ON THE DIURNAL 



sin Au, sin 6u in terms of sin tt d'^Vi / rf/x^. This gives, by substitution and integration 
with respect to X, 

^ Y(3)/a =: (^ 0-0401 Tg^ + 0*1426 T^^ + 0-2523 T/) cos 3X 

+ (0-0319 18^-0-093.6 T/ - 0*5577 T^^) sin 3X 

+ (0*00033 Tg^ + 0-00279T/ -- 0-00411 T^-^) cos 4X 

+ (— 0-00075 Tg^+ 0-00078 T/ - 0-01978 Tg*) sin 4X. 

We have obtained thus finally, 38 coefficients in the expansion of V/a, which, for the 
sake of reference, I collect into Table VIL 

In this Table C«^ is the coefficient of T«^ cos mX. 



Q m 



J? 



}i 



Tn^ sin mX. 



Table VII. 

The unit of force is C.G.S. 10~^. 



02^ 

S,i 
S3I 

Sgi 

S71 


48-00 
81-52 

- 6-66 
+ 7-62 
-f 1-09 
4- 5-78 
4- 1-57 

+ 24 09 
+ 35-59 

- 8-45 
+ 2-39 

2-41 
1-64 
0-34 


G,' 

S 3 
^> 

S/ 


- 2-948 

7-936 

0-630 

0-341 

+ 0-200 

+ 0-407 

+ 10-764 
+ 4-715 
0-214 
0-731 
+ 0-293 
+ 0-264 


C4^ 

S s 

*^4 

s/ 

S33 

S,4 

s,^ 


- 0-252 

0-143 

+ 0-040 

+ 0-00411 
0-00279 
0-00033 

-f- 0-578 

+ 0-094 

0-032 

-f 0-01978 

0-00078 

+ 0-00075 



In order to show how far these numbers correctly represent the forces which have 
been made use of in their computation I have calculated backwards from the potential 
the force to geographical West. Table VIII. exhibits the results, showing by comparison 
the coefficients calculated from the formulae for the potential with those obtained 
directly by observation. It will be seen that the agreement is satisfactory, except 
for the coefficients a^ and 63. 

In figs. 1, 2, 3, and 4 the curves for the mean of the year are plotted. In these 
curves only the first four terms a^ cos t -{- b^ sin t -{- a^ cos 2t + b^ sin 2t have been 
used. The numbers from which the curves are plotted are given in Table IX. 

The curves are seen to be almost identical, and having, therefore, obtained an 
expression for the potential which correctly represents the observed West force, we 
may turn to the main object of this inquiry and calculate the vertical forces. 

The calculation of the various components from the potential involves the knowledge 
of the differential coefficients of zonal harmonics and of the tesseral harmonics for the 
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Bombay. 

Comparison between calculated and observed curve of West force. TJbe abscissse denote astronomical 
time, the ordinates the force to geographical West, the nnit of force being C.Gr.S. 10~^. 

Observed curve, white line. Calculated curve, dotted line. 



Fig. 2. 




Lisbon. 

Comparison between calculated and observed curve of West force. The abscissae denote astronomical 
time, the ordinates the force to geographical West, the unit of force being C.Gr.S. 10""^. 

Observed curve, white line. Calculated curve, dotted line. 

MDOCCLXXXIX. — A. 3 Q 
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Pig. 3. 




Geeenwioh. 

iCompaiison between calculated and observed carye of West force. Tbe abscissse denote astronomical 
time, the ordinates the force to geographical West, the nnit of force being OiGr.S. 10-^, 

Observed cnryOj wMte line. Calculated cnrve, dotted line. 



ig, 4. 




St. Petirsbueg. 

Comparison between calcnlated and observed curve of West force. The abscissge denote astronomical 
time, tb# ordinates the force to geographical West, the nnit of force being O.G.S. lO""*"*. 

Observed curve, white line. Calculated curve, dotted line. 
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four observing stations. The numbers which I have computed for this purpose are 
OTven in Table X. 



In that table T/ stands for 



9T-> 



#1 



Q olXX wj 



where (i = cos u and P^ is the igonal harmonic of degree p. 



Table VIIL 



•Comparison between observed and calculated coefficients of West 

force. 



Bombay . . 



Lisbon . . . 



Greenwich . . 



St, Petersburg 



Winter , 
Mean . , 

S,-j YYi YY\ PT* 

Winter . 
Mean . . 
Summer . 
Winter 
Mean . . 
Summer . 
Winter . 
Mean . . 



a,. 



Observed. 



+ 58 

+ 29 
+ 74 
+ 87 
+ 80 
+ 85 
+ 105 
+ 95 
+ 100 



Oalctt- 
lated. 



+ 58 

+ 27 
+ 88 
+ 52 
+ 70 
+ 91 
+ 95 
+ 93 
+ 79 
+ 121 
+ 100 



^1. 



Observed. 



+ 116 

+ 32 
+ 74 
+ 225 

"[* X JLij 

+ 172 

+ 233 
+ 115 
+ 174 
+ 225 
+ 96 
+ 161 



Calcu- 
lated. 



+ 124 

+ 82 
+ 232 
+ 112 
+ 172 

+ 255 
+ 113 
+ 184 
+ 229 

+ 161 



«•>. 



Observed, 



+ 132 


+ 


4 


+ 


68 


+ 108 


+ 


31 


+ 


70 


+ 


83 


+ 


1 


...- 


.42 


+ 


99 


■■■www 


15 


-- 


42 



Calcu- 
lated. 



+ 124 

+ 8 

+ 66 

+ 90 

+ 10 

+ 50 

+ 85 

- 5 
+ 40 
+ 92 

- 6 

+ 4o 



h^ 



Observed. 



Calcu- 
lated. 



+ 


60 


+ 


23 


+ 


41 


+ 158 


+ 104 


+ 131 


+ 135 


+ 


87 


+ in 


+ 115 


+ 


69 


+ 


92 



+ 73 

+ 39 
+ 56 
+ 172 
+ 105 
+ 138 
+ 149 
+ 104 
+ 127 
+ 100 
+ 72 
+ 87 



Bombaj . . . 
Lisbon . . . 
Greenwicb . . 
St. Petersburg . 


Mean of year 

5? 


ag. 


h' 


a4. 


h* 


Observed. 


Calcu- 
lated. 


Observed. 


Calcu- 
lated. 


Observed. 


Calcu- 
lated. 


Observed. 


Calcu- 
lated. 


+ 71 
+ 63 
+ 41 
- 12 


+ 19 

+ 108 

+ 82 

- 7 


+ 8 
+ 57 
+ 44 
- 8 


- 24 
+ 88 
+ 74 
+ 16 


+ 24 
+ 24 
+ 17 
- 18 


+ 20 
+ 21 
+ 31 

- 20 


- 14 

+ 8 
+ 14 

- 1 


- 12 

+ 6 
+ 21 

+ 27 



S Q 
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Table IX. — Comparison between calculated and observed variation of force to 
geographical West. The unit of force is C.G.S. 10"^. The first four terms 
only of the series a^ cos t -{- 6^ sin ^ + % cos 2t + b^ sin 2t have been taken into 
account. 





Bombay. 


Lisbon. 


Greenwich 


St. Petersburg, 


Observed. 


Calculated. 


Observed. 


Calculated. 


Observed. 


Calculated. 


Observed. 


Calculated. 


A stron OTTi 1 cf^l 


















time. 


















- 


ai + 29 


+ 27 


+ 80 


+ 70 


+ 95 


+ 93 


+ 96 


+ 100 




\ + 74 


+ 82 


+ 172 


+ 172 


+ 174 


+ 184 


+ 161 


+ 161 




a, + 68 


+ 66 


+ 70 


+ 50 


+ 42 


+ 40 


+ 42 


+ 43 




h + 41 


+ 56 


+ 131 


+ 138 


+ 111 


+ 127 


+ 92 


+ 87 


Hour. 





















+ 97-0 


+ 93-0 


+ 150-0 


+ 120-0 


+ 137-0 


+ 133-0 


+ 138-0 


+ 143-0 


1 
JL 


+ 126-6 


+ 132-5 


+ 247-9 


+ 224-4 


+ 228-8 


+ 235-6 


+ 216-8 


+ 219-0 


2 


+ 131-6 


+ 145-9 


+ 303-7 


+ 291-1 


+ 286-4 


+ 302-5 


+ 264-3 


+ 263-9 


3 


+ 113-8 


-F 133-1 


+ 309-2 


+ 309-1 


+ 301-2 


+ 322-9 


+ 273-7 


+ 271-5 


4 


+ 80-1 


4- 100-0 


+ 267-4 


+ 278-5 


+ 273-3 


+ 295-8 


+ 246-1 


+ 243-2 


6 


-r 40-6 


+ 57-0 


+ 191-8 


•+ 210-0 


+ 211*8 


+ 230-7 


+ 190-0 


+ 187-7 


6 


+ 6-0 


4- 16-0 


+ 102-0 


+ 122-0 


+ 132-0 


+ 144-0 


+ 119-0 


+ 118-0 


7 


15-4 


- 13-0 


+ 19-4 


+ 35-8 


+ 51-6 


+ 55-5 


+ 48-2 


+ 48-9 


8 


19-9 


- 24-0 


39-4 


30-5 


13-9 


- 17-2 


- 9-3 


- 7-4 


9 


9-2 


- 171 


- 66-0 


- 65-9 


55-2 


- 62-7 


- 46-1 


-^ 43-9 


10 


+ 10-4 


+ 2-1 


- 61-7 


69-1 


-^ 70-4 


78-5 


61-3 


59-9 


11 


+ 29-6 


+ 24-3 


37-7 


- 48-8 


- 65-8 


71-0 


60-6 


- 61-2 


12 


-f 39-0 


■t- 39-0 


10-0 


- 20-0 


-- 53-0 


53-0 


- 54-0 


- 57-0 


13 


+ 32-2 


-f 37-9 


+ 4-3 


+ 0-2 


- 45-0 


39-4 


- 52-0 


- 57-6 


14 


+ 7-4. 


-f 171 


6-9 


- 2-1 


- 52-2 


42-5 


- 62-9 


- 70-3 


15 


- 31-8 


- 21-1 


47-2 


33-1 


-- 79-2 


- 68-9 


- 89-7 


- 97-5 


16 


77-1 


69-0 


- 110-6 


89-5 


- 123-1 


- 115-8 


128-7 


- 135-6 


17 


- 117-4 


115-4 


^ 182-0 


- 158-6 


- 173-6 


- 172-9 


- 170-8 


- 175-1 


18 


- 142-0 


- 148-0 


242-0 


- 222-0 


•- 216-0 


224-0 


- 203-0 


- 204-0 


19 


143-4 


- 157-4 


271-6 


- 260-4 


- 235-4 


- 251-7 


-213-0 


- 210-3 


20 


- 119-1 


- 139-0 


-~ 257-4 


- 258-5 


- 220-3 


-- .242-8 


- 192-1 


- 186-2 


21 


72-8 


- 94-9 


-~ 196-0 


- 210-1 


-166-8 


- 191-3 


- 137-9 


- 130-1 


22 


13-4 


- 331 


95-1 


- 119-9 


- 79-8 


- 101-5 


- 56-1 


- 47-7 


23 


+ 47-2 


+ 34-1 


+ 27-9 


- 2-6 


+ 27-6 


+ 13-2 


+ 41-4 


+ 48-6 
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Table X. — Numerical values for the zonal harmonics, their differential 


coefficients, 






and the functions T/. 








Pi . . . . 


Bombay. 


Lisbon. 


Greenwich. 


St. 


Petersburg. 


+ 0-3239 


+ 0-6255 


+ 


0-7824 


+ 


0-8664 


P, . . . . 


- 0-3426 


+ 0-0868 


+ 


0-4183 


+ 


0-6235 


Ps • • • • 


- 0-4009 


- 0-3265 


+ 


0-0238 


+ 


0-3224 


P, . ■ . • 


+ 0-0297 


— 0-4225 


— 


0-2811 


+ 


0-0206 


Pb • •• • 


+ 0-3380 


- 0-2144 


, — 


0-4149 




0-2258 . 


Pe .... 


+ 0-1760 


+ 0-1061 


— 


0-3620 


— 


0-3754 


dFcildjui . . . 


+ 0-9718 


+ 1-8764 


+ 


2-3473 


+ 


2-5963 


dJ^^jd/i . 


- 0-7] 31 


+ 1-4341 


+ 


3-0914 


+ 


4-1174 


dFJd;^ , . . 


- 1-8346 


- 0-4089 


+ 


2-5142 


+ 


4-8529 


dF^ld/i . . . 


- 0-4457 


- 2-3681 


+ 


0-5618 


+ 


4-3029 


dF^jdjui , , . 


+ 1-8839 


- 2-7679 


-- » 


2-0498 


+ 


2-3692 


dFqldjUL . . . 


+ 1-8421 


- 0-9882 


— 


4-1305 




0-5776 


d^Foldfi^ . . 


+ 3-0000 


+ 3-0000 


+ 


3-0000 


+ 


3-0000 


d%/d^^ . . 


+ 4-8588 


+ 9-3821 


+ 


11-7364 


+ 


12-9816 


d^'PJdfj? . . 


- 1-9916 


+ 13-0386 


~f~ 


24-6401 


+ 


31-8221 


d^FJdjLL^ . . 


- 11-6529 


+ 5-7017 


+ 


34-3644 




56-6573 


d^FJdjii^ . . 


- 6-8945 


- 13-0104 




30-8207 


+ 


79-1540 


d^Vrildfl^ . . 


+ 12-8380 


-- 30-2804 


+ 


7*nm 


+ 


87-4567 


d^FJdjLL^ . . 


+ 20-7369 


-- 27-8334 


— 


31-1374 


+ 


70-4907 


^PJ^^/^S . , 


+ 34-011 


+ 65-674 


__ 


82-156 


+ 


90-872 


d^FJd/i^ . . 


- 94-170 


+ 128-393 


+ 


460-164 


+ 


714-102 


tZ3Pg/c?;t3 . , 


+ 98-399 


- 325-813 


+ 


575-912 


+ 


2025-953 


^Ps/tZ/^s . . 


+ 306-10 


+ 591-07 


+ 


739-39 


^ mm 


817-84, 


d%id/ii^ . . 


-918-11 


+ 2260-17 


+ 


6721-52 


+ 10101-17 


d^F^ldjii^ . . 


+ 754-68 


-3278-66 


+ 16512-03 


+ 44542-37 


T,i . . . . 


+ 0-9194 


+ 1-4640 


-1- 


1-4620 


+ 


1-3008 


^^^ • • • • 


0-6747 


+ 1-1190 




1-9253 


+ 


2-0629 


T3 


- 1-7358 


- 0-3190 


+ 


1-5658 


+ 


2-4313 


^^ • • • • 


- 0-4217 


- 1-8477 




0-3498 


+ 


2-1558 


l< • • • • 


+ 1-7825 


- 2-1596 


— _ 


1-2766 


+ 


1-1870 




+ 1-7429 


- 0-7710 


— 


2-5725 




0-2894 


'iV .... 


-t 2-685 


+ 1-826 




1-163 


+ 


0-763 


T32 . . . . 


+ 4-349 


+ 5-712 


+ 


4-551 


+ 


3-258 


T/ ..... 


- 1-783 


+ 7-938 




9-556 


+ 


7-987 


±5 . . . . 


- 10-430 


+ 3-471 


H- 


13-327 


+ 


14-221 


Xg . . . . 


- 6-171 


- 7-921 


+ 


11-953 


+ 


19-868 


T^^ . . . . 


+ 11-491 


- 18-435 




2-993 


+ 


21-962 


T32 . . . . 


+ 18-561 


- 16-945 


— 


12-075 


+ 


17-693 


3 


+ 28-801 


+ 31-195 


+ 


19-841 


+ 


11-427 


rp 3 


- 79-745 


+ 60-985 


+ 


111-130 


+ 


89-801 


W 3 

Xg . . , . 


+ 83-326 


- 154-760 


+ 


139-090 


+ 


254-770 


rp 4 

X5 . . . . 


+ 245-23 


+ 219'06 


+ 


111-20 


+ 


61-525 


T^^ , . . . 


--7^5-55 


+ 837-66 


+ 


1010-90 


+ 


636-390 


T,* .... 


+604-62 


-1216-10 


+ 


2483-40 


+ S 


2806-300 
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IV. Comparison of Observed and Calculated Vertical Forces, 

The complete potential can be written down from the value at its surface in the 
usual manner, either on the supposition that the potential is zero at an infinite 
distance, or that it is zero at the centre of the Earth, the first supposition 
corresponding to the hypothesis that the seat of the magnetic variation is outside 
the Earth. If Y^ is a surface harmonic of degree n occurring in the expansion of 
V/a, the solid harmonic will either be YJ^rjaf^ or Yn{rlaY^^^'^^\ The vertical force 
is given by dY/dr^ as the force is considered positive when it acts downwards. At 
the surface, therefore, we have a term for the vertical force which is either nY^ or 
^{n + l)Y,. 

Before we proceed to discuss the comparison between the observed and calculated 
values of the vertical force, a few words are necessary regarding the available 
observations. 

The only station for which we have complete records for 1870 is Lisbon. It is 
therefore impossible to obtain a satisfactory series for the vertical force which 
would give us, if our information was more complete, directly the two terms, one 
due to outside, the other due to inside, effects. But I shall show that even from the 
existing data we can draw important conclusions. 

At Bombay no vertical force determinations are published, as far as I know, before 
1873, when the magnetograph came into operation; but we have complete records 
between 1873 and 1878. During these years the general type of the vertical force 
remained practically the same, only the range varying. Figs. 5 and 6, for instance, 

Fig. 5. Eig. 6. 
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are tracings of the curves given by Mr. Chambers in the years 1873 and 1877, the 
two years which differ most in range. The other years show curves varying between 
these two extremes. As far as the general form of the curve is concerned, we cannot 
go far wrong, therefore, if we make use of the 1873 observations, especially as the 
horizontal components show no marked difference (except as regards range) in 1873 
and 1870. 

Similar remarks apply to Greenwich. Although published records exist for 1870, 
there is an uncertainty about the temperature correction, which makes the vertical 
force observations previous to 1882 useless for our present purpose. 

Table XI. will give an idea of the changes which the vertical force variation is 
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subject to in various years. The numbers are copied from the published records of 
the Greenwich Observatory. 



Table XI. 
Coefficients in the series 

V^ = Ci sin {t + a) + Cg sin {2t + ^) + c^ sin (3^ + y) + ^4 ^^^ (^^ + y) 

where changes in V^; represent changes in the vertical force, the unit being 
•00001 of the whole vertical force ; t being the time from midnight. 





■ 


o 


Cg. 


^. 


C3. 


7- 


^4- 


a. 






o 












1883 


14-3 


148-13 


13-1 


266-58 


5-3 


89-60 


1-3 


293-20 


1884 


14-8 


139-33 


11-7 


272-00 


5-5 


95-52 


21 


289-49 


1885 


130 


137-50 


11-7 


265-35 


5-1 


83-28 


1-5 


281-04 


1886 

■ 


12-6 


160-58 


11-9 


268-38 


4-0 


94-22 


1-2 


297-50 



The principal discrepancy here occurs in the angle a, which would give a difference 
of phase between 1885 and 1886 of about 23°, or about an hour and a half; the 
angle a differs much at Greenwich during different months ; the phases of the other 
terms show practically no difference. In comparing the observed and calculated 
values, I have taken the year 1884, as during that year the range of the declination 
needle was greatest, and corresponded most nearly to that of 1870. 

At St. Petersburg, also, the results for 1870 are not corrected for temperature, and 
there is reason to beheve that this has affected the observed type cjonsiderably. I 
have therefore compared the results of calculation with observation, also, for the year 
1878, in which year the temperature corrections have been taken into account. 

We must remember, then, in the comparison between the observed and calculated 
values, that the greatest weight must attach to the Lisbon observations, and that 
comparatively little value is to be given to St. Petersburg, as I have no records at my 
disposal which would enable me to judge how far the type of vertical force ^ there 
varies from year to year. 

Tables XII. and XIII. give the values of the coefficients a^, b^, a^, h^, &a, as calculated 
on the hypothesis that the disturbing force comes from the inside of the Earth. 
Tables XIV. and XV. give their values calculated on the hypothesis that the 
disturbing force is outside. In both cases the observed numbers are given for 
comparison. I have calculated the values for Lisbon separately for the winter and 
for the summer months ; in the other cases the values for the mean of the year only 
have been taken. The years attached to the observing stations refer to the date of 
the observations ; the calculated values all belong to 1870. 
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Table XII. — Comparison between the observed and calculated coefficients for the 
vertical force, a^, &i, %/^2 ^^ the hypothesis that the disturbing force is inside 
the Earth. 



Bombay . . . . . . 

Lisbon, summer . . . 
„ winter .... 

Greenwich, 1884 . . . 
St. Petersburg, 1870 . . 

1878. . 


%. 


K 


1 
^2- 


5o. 


Calculated. 


Observed. 


Calculated. 


Observed. 


Calculated. 


Observed. 


Calculated. 


Observed. 


+ 219 
+ 580 
+ 336 
+ 350 
+ 177 


- 42 
-196 

- 135 

- 42 
+ 114 

8 


- 57 
~~206 
-149 
-190 
-154 
* . 


+ 9 
+ 66 
+ 55 
+ 49 
+ 125 
+ 29 

■ 


4- 79 
+ 381 

+ 277 
+ 139 

+ 62 

. « 


■ 

- 21 

-201 

-103 

51 

- 59 
24 


-152 

- 87 

- 19 

- 34 

- 45 

• • 


+ 28 
+ 26 
+ 1 

+ 2 
38 

- 2 



Table XIII. — Comparison between the observed and calculated coefficients ag, 63, a^, 64 
for the vertical force on the hypothesis that the disturbing force is inside the 
Earth. 



Bombay . . . 
Lisbon . . . 
Greenwich. . . 
St. Petersburg . 



■ 


Jg. 


a^. 


54. 

i 


Calculated. 


Observed. 


Calculated. 


Observed. 


Calculated. 


Observed. 


Calculated. 


Observed. 


- 73-4 
+ 156-1 
+ 85-5 
+ 12-0 


- 10 

- 81 
.- 24 

- 12 


- 42-0 
-171-4 

- 88-0 

- 17-4 


+ 35 
+ 12 

+ 3 

- 8 


- 20-4 
+ 9-2 
+ 28-1 
+ 22-3 


+ 3 

- 21 

- 8 

- 18 


- 29-0 

- 221 

- 38-2 

- 33-3 


+ 16 
+ 11 
+ 3 
- 1 



Table XIV.— Comparison between the observed and calculated coefficients a^, h^, a^, h^ 
on the assumption that the disturbing force is outside the Earth. 



Bombay . . . . • • 
Lisbon, summer ... 
„ winter .... 
Greenwich, 1884 . . . 
Sfc. Petersburg, 1870 . . 

1878, . 


«i. 


h. 


%. 


^2* 


Calculated. 


Observed. 


Calculated. 


Observed. 


Calculated. 


Observed. 


Calculated. 


Observed. 


-141 

398 

- 248 
235 

- 97 

• • 


42 
196 
135 

- 42 
+ 114 

- 8 


+ 31 
+ 137 

+ 108 
+ 132 
+ 104 

• • 


+ 9 
+ 66 
+ 55 
+ 49 
+ 125 
+ 29 


53 
315 
235 
110 

- 40 

• ♦ 


- 21 

-201 

103 

51 

- 59 

- 24 


+ 121 
+ 52 
+ 15 
+ 21 
+ 35 
» • 


+ 28 
+ 26 

+ 1 
+ 2 
- 38 

2 
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Table XV. — Comparison between the observed and calculated coefficients ag, 63, %, b^ 
on the assumption that the disturbing force is outside the Earth. 



Bombay . . . 
Lisbon . . 
Greenwich . 

Sfc. Petersburg . 



%. 



Calculate<3. 



+ 66 
-134 
- 70 

H 



Observed. 



10 
81 
24 
12 



h,. 



Calculated. 



4- 125 

+■ 38 
+ 16 



Observed. 



+ 35 

+ 12 

-f 3 

- 8 



^4. 



Calculated. 



+ 18 

- 8 

- 25 

- 20 



Observed. 



+ 3 

- 21 

- 8 
~ 18 



h. 



Calculated. 



+ 24 

-f- 18 

-h 33 

+ 27 



Observed. 



+ 16 

+ 11 

_|_ ^ 

- 1 



Confining ourselves in the first place to the first four coefficients, we find that out 
ot twenty coefficients eighteen have the wrong sign on the hypothesis of an internal 
cause, while only two have the w^rong sign on the hypothesis of an external cause, 
and those two belong to St. Petersburg, to which station, as was pointed out, we cannot 
attach much value. If instead of the numbers given for 1870 we take the numbers 
given at the same station for 1878, the agreement becomes better, even for St. 
Petersburg. The coefficients a^, 63, %, 64 are of course more uncertain ; but even 
here the evidence is strongly in favour of the external cause. Out of twenty 
coefficients seventeen agree in sign with that hypothesis. 

A better comparison can perhaps be obtained in a different way : the two terms 



can be written 



a,i cos nt + bfi sin nt 
Tn cos n{t -— tf), 



where r^ is the amplitude of the oscillation, and t^ the time at which the maximum 
elongations take place. 

Tables XVI. and XVII. contain the results for tn, and from these tables I think it 
will clearly appear that the phase of the vertical force completely agrees with the 
assumption of an external cause and completely disagrees with the assumption of an 
internal cause. 

For Lisbon, our principal station of comparison, the phase in Table XVI. agrees 
for both the diurnal and semi-diurnal variation within four minutes of time. For 
Bombay the diurnal variation agrees within three minutes, and the semi-diurnal 
variation within 36 minutes. For Greenwich, the semi-diurnal variation, which we 
have seen differs little from year to year, agrees closely, while the diurnal variation 
shows a greater difference. In all these cases the phases, as calculated in Table XV., 
are in as great a disagreement as possible. St, Petersburg gives less decisive results, 
but they still go in the same direction, especially if we take the observations of 1878 
to represent the type of vertical variation. 
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The amplitudes of the variation are given in Table XVIII. where the second 
column gives the calculated values on the hypothesis which we must reject, while 
the third column gives the same numbers calculated on the assumption of an 
external cause. 

The calculated numbers are much larger than the observed ones, which is a fact 
requiring explanation. But, the range of vertical force differing from year to year, 
we must confine ourselves in our reasoning principally to the Lisbon observations. 

The changes in the range of different years, as far as we have any observations, 
are not so considerable, however, as to account for the discrepancies at Bombay or 
Greenwich, and we must conclude, I believe, that also for those stations a considerably 
smaller range is obtained in the observed than in the calculated forces. 



Table XVI. — Observed and calculated values of the coefficients t^ and t^ of vertical 
force, when expressed in the form r^ cos {t — t^) + r^ ^^^ 2 {t — t^) on the sup- 
position that the disturbing force is inside the Earth. 



Bombay 

Lisbon 




t,. 




t,. 


Calculated. 


Observed. 


Difference. 


Calculated. 


Observed. 


Difference. 


h. m. 

23 02 
22 35 
22 06 
21 16 

• • 


h. m. 

11 13 

10 4.0 
8 42 
3 10 

7 05 


h. m. 
+ 11 49 

+ 11 58 
11 57 

- 5 54 
9 49 


h. m. 

9 55 

11 42 

11 32 

VlO 48 

• • 


h. m. 

4 23 

5 50 

5 56 
7 05 

6 12 


h. m. 

+ 5 32 . 
+ 5 52 
+ 5 36 
+ •3 43 
+ 4 36 


Greenwich . 
St. Petersburg, 

55 


• * • • 

1870. . 
1878 . . 



Table XVII. — Observed and calculated values of the coefficients t^ and t^ when 
expressed in the form 7\ cos {t — t^) + r^ cos 2 (^ — t^) on the supposition that the 
disturbing force is outside the Earth. 



>■■■■ "•■" — -,,■■,1.1--...— .,,-. -■■-, -— ■— ' —.....-..- 

Bombay 

Lisbon 


. ^i- 




■ ^2. 




Calculated. 


Observed, 


Difference . 


Calculated. 


Observed. 


Difference. 


h. m. 
11 10 
10 37 

10 03 

8 52 

* • 


h. m. 
11 13 

10 40 
8 42 
3 10 

7 05 


h. m. 

-0 03 

03 

+ 1 21 

+ 5 42 
-^1 47 


h. m. 

3 47 
5 46 
5 38 

4 38 
• . 


h. m. 

4 23 

5 60 

5 56 
7 05 

6 12 


h. m. 
-0 36 
-0 04 
-0 18 
-2 27 
-1 34 


Green wicli , 
St. Petersburg, 


• • • • 

1870 . . 
1878 . . 
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Table XVIII. — Observed and calculated values of r^ and r^ in the expression 

r-y cos [t — t-^ + ^2 <^os 2 [t — t^) for vertical force. 



Bombay 

Tjisbon ...... 

Grreenwicli 

St. Petersburg, 1870 . . 

1878. . 


n. 


^2. 


Calculated 
from inside. 


Calculated 
from outside. 


Obseryed. 


Calculated 
from inside. 


Calculated 
from outside. 


Obseryed. 


226 
491 

398 
235 

• • 


144 
346 

269 
142 

• • 


43 
176 

65 
169 

30 

■ 


171 
333 
143 

77 

■ * * . 


132 

277 

112 

53 

• • 


35 

153 

51 

71 

24 



The agreement between the calculated and observed curves for vertical force is best 
seen from the graphical representation given in figs. 7, 8, 9; 10. In the curves, the 
diurnal and semidiurnal variation only have been taken into account. It is seen 
how at Bombay, Lisbon/ and Greenwich the observed curves are almost identical in 
shape with the curves calculated on the hypothesis of an external cause, if the range 
is reduced in a proper ratio. The disagreement with the curves calculated on the 
other alternative is complete, a maximum occurring where a minimum should occur, 

'■■■■■ ■ ■ -^ ' - 

and vice versd. At St. Petersburg, although the agreement is much worse, the 
1878 curve, which has been corrected for temperature, follows pretty closely the shape 
of the calculated curve. The numbers from which the curves have been drawn are 
given in Table XIX. 



Fig, 7. 




Bombay. 
CoDaparison between calculated and observed curve of vertical force. The abscissa denote astronomical 

time, tbe ordinates vertical force, the unit of force being O.G.S. 10"^. 
Observed curve, white line. Curve calculated on hypothesis of outside force, dotted line. Curve 

calculated on hypothesis of inside force, white and black line. 
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Fig. 8. 
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Lisbon. 



Comparison between calculated and observed curve of vertical force. The abscissae denote astronomical 

time, the ordinates vertical force, the unit of force being C.G.S. 10~^. 
Observed curve, wbite line. Curve calculated on hypothesis of outside force, dotted line. Curve 

calculated on bypotbesis of inside force, white and black line. 
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Greenwich. 

Comparison between calculated and observed curve of vertical force. The abscissae denote astronomical 

time, tbe ordinates vertical force, the unit of force being O.G.S. 10"^. 
Observed cnrve, white line. Cnrve calculated on hypothesis of outside force, dotted line. Curve 

calculated on hypothesis of inside force, white and black line. 
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Comparison between calculated and observed curve of vertical force. The abscissae denote astronomical 

time, the ordinates vertical force, the unit of force being C.G.S. 10~^. 

Observed curves, white line — I. for 1870, II. for 1878. Curve calculated on hypothesis of outside force, 

dotted line. Curve calculated on hypothesis of inside force, white and black line. 
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If we, then, take it as proved that the primary cause of the variation comes to us 
from outside the Earth's surface, we are led to consider that a varying magnetic 
potential must cause induced currents within the Earth, if that body is a sufficiently 
good conductor. These induced currents might be the cause of the apparent reduction 
in amplitude. As my colleague, Professor Lamb, had given considerable attention to 
the problem of cuiTents inside a conducting sphere, I consulted him, and he gave me 
the formulae by means of which the induced currents can be calculated. His investi- 
gation is added as an Appendix to this paper. 

T. Discussion of effects due to Currents Induced in the Inside of the Earth. 

I shall assume, then, for the present, that there is a periodic magnetic disturbance 
having its cause outside the Earth, and being probably due to electric currents in our 
atmosphere. Currents will be induced within the Earth, and we must now discuss 
what the effect of these currents will be, and whether they will account for the reduc- 
tion in amplitude of the vertical forces which the observations show. 

The varying potential can be expressed as a sum of terms of the form 

O^ cos {pt 4" ^)j 

where H is a solid harmonic of degree n. Professor Lamb's formulae allow us to 
calculate for .each value of % and for each value of p, the magnetic effect due to the 
induced currents, on the supposition that the specific resistance p of the Earth is 
uniform. The forces due to these currents will have a different amplitude and a 
different phase from the original forces, and it is the resultant effect which we observe 
in the diurnal variations. The general effect will be to increase the horizontal com- 
ponents and to diminish the vertical component. The difference of phase will be the 
same for all components, provided we give a different sign to the amplitude of the 
vertical components of the inside and outside currents respectively. Otherwise the 
difference of phase of the vertical component will be greater by two right angles than 
the difference of phase of the horizontal components. 

If one of the horizontal forces and the vertical force due to the solid harmonic of 
positive degree n are written 

a cos {pt + X) and h cos {pt + X -f e), 

the corresponding components due to the induced potential of negative degree will be 
of the form 

c'a cos [pt +^ + ^) and ch cos {pt + X + € + ^)- 

Table XX. gives the coefficients d, c, and a for given values of the specific resis- 
tance p, if ^ = 2. The value of S has the same meaning as in Professor Lamb's 
paper, and is connected with p by the equation 

pS — AirpW 
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where p is equal to 27rm/T, and m is equal to 1 for the diornal variation, and equal to 
2 for the semidiurnal variation. 

Table XXI. gives the same quantities for n = 4. 

An example may render the use of the Table more inteUigible, Let the Earth, for 
instance, have a uniform specific resistance, v^hich in C.G.S. units is 1*23 X 10^^, and 
consider that term of the magnetic potential which, on the surface of the Earth, has 

the form a • ^^^ 

A sm u cos u cos -^ ? 

X 

where T is the time of the revolution of the Earth. This term alone represents 
fairly well the characteristic features of the diurnal variation. Here, as ^ = 2, we 
may use Table XX. The value of S coiTesponding to the assumed resistance is 30? 
and we find c = — *5 approximately, which means that the vertical force due to the 
induced currents has half the amplitude of the vertical force due to the primary 
currents ; also, the difference in phase is 41° ; as the sign of the amplitude is changed, 
this means that the minimum of the vertical force due to induced currents takes place 
not quite three hours after the maximum of the corresponding primary force. Similarly 
the horizontal force due to the inside current is only one-third of the corresponding 
horizontal force due to the outside potential, and here the maximum due to the 
secondary currents takes place nearly three hours after the maximum due to the 
secondary currents. To get similar numbers for the semidiurnal variation we should 
have to put m = 2, and find c, c\ and <x for S = 60, because p = 1*23 X 10^^ 
= 2 X 6'15 X 10^^, and looking up 6 '15 X 10^^ in the last column we should find the 
corresponding number in the first column to be 60. 

Table XX. — Comparison of the magnetic forces due to a system of varying electric 
currents outside the Earth, with the forces due to the currents induced inside 
the Earth, n •=■ 2, 



I 


c. 

Eatio of normal 

forces due to 

secondary and primary 

variation. 


Ratio of tangential 

farces due "to 

secondary and primary 

variation. 


a. 

Difference of phase.. 


Corresponding value of p. 


1 


- -02854 


+ *-01903 


87 27| 


3-70 X 1014 X m. 


2 


'05688 


•03792 


84 55| ^ 


1-85 X 101* 


3 


•08484 


•05656 


82 25 


1-23 X 1014 


4 


•11226 


•07484 


79 57 


9-25 X 1013 


6 


•13895 


•09263 


77 32| 


7^40 X 1013 


6 


•16479 


•10986 


75 10| 


6-17 X 1013 


7 


•18968 


•12645 


72 53 


529 X 1013 


8 


•21B54 


•14236 


70 41 


4-63 X 1013 


9 


•23633 


•15755 


68 32 


4-11 X 1013 


10 


•25799 


•17199 


66 29 


3-70 X 1013 


20 


•41771 


•27847 


• 50 32 


1-85 X 1013 


30 


•50520 


•33680 


41 03 


1^23 X 1013 


40 


•56158 


•37439 


35 11 


9-25 X 1013 


60 


•59864 


•39909 


31 11 


7-40 X 1012 


100 


•69949 


•46632 


21 33 


3^70 X 1012 
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Table XXI.-— Comparison of the magnetic force due to a varying potential repre- 
sented by a solid harmonic of degree 4 with the corresponding forces due to 
currents induced inside the Earth. 



d. 


c. 

Eatio of normal 

forces due to 

secondary and primary 

variation. 


Ratio of tangential 

forces due to 

secondary and primary 

variation. 

■ 


a. 

Difference of phase. 


Corresponding value of p. 


1 

10 

100 


- -0101 
•0992 
•5149 


+ -0081 
•0794 
•4119 


89 01 
80 21 

37 50 


3-70 X 10^4 X m. 
3-70 X 1013 
3-70 X 10^3 



Tables XX. and XXI. cannot, however, be used to compare our calculated and 
observed results, but form only an intermediate step. 

We observe on the Earth the resultant of the outside and inside effect, and we 
have calculated the vertical force on the assumption that the whole horizontal force 
is due to outside effect. 

In fig. 11, let OHj^ represent that part of the horizontal force which is due to the 

Fig. 11. 




outside effect; and OH^ the corresponding force in phase and amplitude which is due 
to the induced effect. The observed horizontal force will be the resultant OK. Let 
OVi be the magnitude of the vertical force due to the outside currents, and OV^ the 
vertical force in phase and amplitude due to the induced effect. The observed 
vertical force will be ON. If w^e calculate the vertical force on the assumption that 
the resultant OK is due entirely to the outside, we should obtain a force OE, such 
that OR : OK : : OV^ : OHi. Our calculated value of vertical force will be OR, and 
our observed value ON. We require the ratio of lengths of these lines and the angle 
between them. 

MDCCOLXXXIX. — A. 3 S 
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In the triangle ONE. we know, from Table XX., the ratios 



V^R : OVi = HiK : OH, = OH3 : OH, = c, 

V-^N : OVi = OV3 : OV, = c, 



also the angle 



Write 



NViO = V1OH3 = a. 



NOVi = 6 • ROYi = ^' ; ONVi = 4, • ORVi = f . 



The triangles OV^N and OV^R give us the equations 



tan I (</> - ^) = --™"^ cot ^ 



+ c 



2 



(^ ^-' ^) 



77 



06. 



^^' -|- ^^^ zz=L a. 

These equations determine and ^', and hence the required angle y = ^ + ^'' ; also 
OE : ON = sin (^ : sin ^'. 

In Table XXII. the angle y and the ratio ON : OE : = r have been calculated 
for n = 2. Table XXIII. gives the corresponding quantities if the inducing solid 
harmonic is of degree 4. 



Table XXII. —Comparison between resultant vertical force as regards magnitude and 
phase when induced currents are taken into account and vertical force calculated 
on the assumption that the whole is due to an outside effect. The inducing 
potential is a solid harmonic of degree 2. 



a. 


Reduction in 
amplitude. 


Change of pliasc. 
7- 


P' 


1 


•9981 


o / 

2 43 


3-70 X 101^' 


5 


•9565 


13 02 


7-40 X IQi^^ 


10 


•8589 


23 10 


3-70 X 1013 


20 


•6705 


34 03 


1-85 X 10^3 


30 


•5516 


38 12 


1-23 X 1013 


40 


•4762 


40 16 


9-25 X 1013 


50 


•4261 


41 10 


7^40 X IQi^ 


100 


•3004 


43 08 


3-70 X 1012 
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Table XXIIl. — Comparison between resultant vertical force, as regards magnitude 
and direction when induced currents are taken into account, and vertical force 
calcnlated on the assumption that the whole is due to an outside effect. The 
inducing potential is a solid harmonic of degree 4. 



a. 


Reduction in 
amplitude. 


Change of phase. 


P- 


1 

10 

100 


•9997 
•9723 

•4982 


O I 

1 06 
10 06 
38 49 


3-70 X 1014 
3-70 X 1013 
3-70 X 1013 



The observed amplitude of the vertical force at Lisbon is about one half of its cal- 
culated value. If the conductivity of the Earth was such as to produce this reduction 
in amplitude, it is seen from Tables XXII. and XXIII. that the phase would be altered 
about 40°, while in reality there is a remarkable agreement in phase. If the con- 
ductivity is so small as to leave the resultant phase practically unaltered, as observation 
tends to show, the amplitude also should not be sensibly altered. There is, therefore, 
no uniform conductivity of the Earth which can make the observations agree with 
the calculation. Such an agreement, however, can be easily brought about, as 
Professor Lamb has suggested to me, if the conductivity of the inside of the Earth is 
larger than the conductivity of the upper layers. It is extremely probable that this 
is really the case. The bulk of the outside layer of the Earth, except in so far as it is 
water, is made up of material which in its ordinary condition is non-conducting ; but we 
know that some of the siHcates begin to conduct at temperatures above 200"^ C, and, 
generally speaking, insulators lose their insulating powers at high temperatures. With- 
out regard even to the quantities of metallic matter that may be stored inside the Earth, 
there is nothing improbable in the supposition that its conductivity increases towards 
the inside. If the bulk of the observed induced effect is due to currents in a fairly 
conducting inner sphere, the calculated phase would be that due to good conducting 
matter, and would not differ from the observed value, while the reduction in amplitude 
might yet be sufficient to account for the observed facts. In order to give a better 
idea of the kind of conductivity which is required to produce a certain change of 
phase, it may be stated that for the purest distilled water obtained by Kohlrausch 
p would be about 1-4 X 10^^ Such water is, as is well known, a very bad conductor, 
and, according to our Tables, if the whole Earth was made up of matter which conducts 
as badly, there would be no currents in the Earth induced by the diurnal variation of 
sufficient intensity to affect our magnetic needle sensibly. Ordinary rain water, how- 
ever, has a specific resistance of about 6 X 10^^ A conducting sphere of the same 
resistance would already produce a retardation in phase of about an hour for the 
diurnal variation if the solid harmonic is of degree 2. For salt water the resistance 

3 s 2 
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may get as low as 4 X 10^. A whole sphere made up of such water would very 
considerably reduce the amplitude of the observed vertical force, and alter the 
resultant phase by 45° nearly. The average conductivity of the Earth, as seen from 
these examples, must be small, although it may be considerable over hmited areas. 
Such hmited areas would principally affect the harmonic terms of higher degrees, and 
we should not consequently expect for them such a good agreement between theory 
and observation. Table XXII. shows tlmt as the resistance of the sphere diminishes 
the retardation of the resultant phase seems to approach a constant value of 45°. 
This can be proved to be quite generally the case. It follows directly from a formula 
given by Professor Lamb in the Appendix, and can also be seen as follows : — When 
the conductivity is good the angle between OHg and OH^ in fig. 11 will steadily 
diminish, and ultimately vanish. OK will ultimately coincide with OH^, and OR 
with OVj ; but the angle between OV3 and OV^ will increase towards 180^, and the 
sides will tend towards equality. Two very nearly equal and very nearly opposite 
forces may have a resultant which is inclined by a finite angle to the forces. To find 
the angle in the limit between OV^ and OV2 we deduce, in the first place, an expres- 
sion for the ratio of the vertical force due to the outside effect to the vertical force 
due to induced currents for good conductivities. 

This ratio depends, as shown by Professor Lamb, on the function 



and 



;^,e (kR) where k^ = , 

X40 = 3•5•••(2'^+1)(^J -^ 



If n is odd, 



/ d Y . . y cosli^ n,n + 1 sinh^ 7i — l,n.n + l,n + 2 cosli^ 

\Ydu ^^ ^ ^^ yn + i "~" 2 hi+^ "' 2.4 t''*'^'^ 

{n —2) . . . (7^ + 3) sinh ^ 

If n is even, we must interchange sinh ^ and cosh ^. 

If ^ is larger, so that e'^ can be neglected, compared to e^, cosh { = sinh C = e^, and 



2n + 1.(2^ + 3) %(;^-l)(/^E) 

which is the ratio of the vertical forces due to the induced and inducing potential 

will become 

n -{■ l,n -\- 2 y_. , n.oi + l,n + 2.n + 3 ^_^ ^ 
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GenerallY in the abo?e expansion we may put therefore 



^•"f = §-*# { COS "-" — % sin "J , J 



so that the ratio of the vertical forces will be 



where 



XjhJY 
X' + iY 



f i> 



2 4 '2.4 4 

(?i. — 1) . • , (fi. + 4) 3^ ^. 

9 A A '^^^ A 

^' ^ # + 1m + 2 ♦ ^ ^_^| ti,# + 1.^ + 2m + 3 . 2w ^^^ 

X — — " ^ ''' , sm . o * *""" -— ■- -_-—— ~| --. sm J o 



'I 



(^ •-. 1) , . /^ 4. 4) , 37r K . 

JU i— _- £,— -Ji ™_4 o-iri — ^ v«-f 

2 4 6 4 * 



X^ = 1 — - ----- COS 7 o""^ + ~ — — — — ~ — ~ — ■ cos — ^ — . . • , 

5f — l.# . ^ ff — 2»w — l.?*.# + 1 .- 2ir *_. , 
Y ■= '~^-— ■■ sm 7 S""^ — -— — — -"-7 — — sm -- o • -^ + . . • 

2 4 2.4 4 

The ratio of the amplitudes^ if S is large^ becomes 

-^^ilr- + ^ I = 1 •-» /2fi 4- 1) (23^-^ 

and the angle between the two vertical forces 

tan""^ :^ — tan""^ w = (2s +' 1) (2S)""^ . 

The resultant of these two forces^ resolved in a direction parallel to either of them, 
is therefore equal to the component which is at right angles, and the resultant will 
consequently be at an inclination of 45*^. 
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The result will not hold, of course, if conduction in the Earth's crast takes place 
chiefly at some distance below the surface, as in that case the vertical force due to the 
induced currents will not tend to become equal to the vertical force due to the 
primary variation. For a bad conductivity we shall always have a resultant vertical 
force sensibly equal to the primary force. If the conductivity increases, the resultant 
will have a different phase from the primary variation, tending towards a difference 
of 45^, if the conductivity is uniform. If the conductivity is not uniform, a maximum 
difference of phase will be reached, which, if the conductivity is still further increased, 
diminishes indefinitely. 

VI. The Magnetic Potential on the Surface of the Earth. 

As it seemed interesting to trace the equipotential lines on the surface of the Earth 
as far as they depend on the diurnal variation, I have calculated the potential from 
the equations [A] and [B]. 

It is necessary, for this purpose, to compute the tesseral harmonics for definite 
points on the Earth's surface. It would seem most natural to choose these points, so 
that they lie on equidistant circles of latitudes, but as tables'''" exist for the zonal 
harmonics in terms of the cosine of the colatitude, I have selected values of these 
cosines so that the corresponding angles should differ as nearly as possible by 10^ 

The values of u and /x = cos u, for which the potential is computed on the Northern 
hemisphere, are given in Table XXIV., u being the colatitude. 

Table XXIV. 

^=:cosu= -98 '94 -87 77 -64 -50 -34 -17 '00 

u^ 11'' 29' 19' 57' 29° 32' 39'' 39' 50^12' 60^00' 70° 07' 80° 13' 90° 00' 

Symmetrical circles of latitude were taken on the Southern hemisphere. 

Taking from Mr. Glaisher's Table the values of P,- corresponding to each of the 
above values of /x,, we obtain the differential coefficients of zonal harmonics by a 
successive application of the formula 

I - ^ = (2i - 1) P,.,, 

The first and second differential coefficients thus calculated are given in Tables 
XXV. and XXVI. 



* 'Report of tlie British Association ' (Sheffield, 1879). 
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Introducing tbese quantities in the equations for the potential, and taking proper 
account of the change of sign of /x in the Southern hemisphere, I have obtained 
Table XXVII., in which the potentials are given for 24 equidistant meridian circles. 
In order to reproduce the daily variations, we must imagine the whole system of 
equipotential lines to revolve round the Earth from East to West ; the time for which 
the potential is given is mean noon for the zero, meridian. It will be remembered 
that the equations for the potential have been derived from the mean summer values 
in the Northern, and mean winter values in the Southern hemisphere. If we want 
to get a symmetrical potential in both hemispheres, we must take the average varia- 
tion for the whole year, or, what comes to the same thing, we may in Table XXVII. 
write down the mean values for two corresponding circles of latitude, one in each 
hemisphere. This has been done in Table XXVII L, where the values are only given 
for the Northern hemisphere. The mean equipotential lines for the year are drawn 
in fig. 12. If we imagine the variable part of the magnetic force to be produced by 
a system of surface currents in a conducting sphere concentric with the Earth, and 
surrounding it, we may, if the potential is known, calculate the distribution of the 
lines of flow. 

If the magnetic surface potential is of the form O,^, when O,^ is a harmonic of degree n, 
the current function ^,, is given by 

n + 1 

SO that the lines of flow are the same as the equipotential lines. This is no longer 
true when the magnetic potential is made up of a number of terms corresponding to 
harmonics of different degrees, for the factor {2n + l)/(^^ +1) will vary for different 
terms, and the resultant current function will therefore no longer be proportional to 
the resultant magnetic potential. 

In our own case, taking the mean values for the whole year, the series begins with 
O2, and the factor {2n + l)/(n + 1) will vary, therefore, only between 5/3 aiad 2. We 
may then, as an approximation, still take the equipotential lines to give us the general 
form of the lines of flow. We conclude that we may imagine the daily variation 
of the Earth's magnetic force to be produced by a system of electric currents in a 
sphere surrounding the Earth, in which the lines of flow are roughly represented in 
fig. 12, the direction being such that at longitude eo"" East the flow is away from the 
equator. 



*' Maxwell, ' Electricity and Magnetism/ vol. 2, p. 281. 
MDCCCLXXXIX. — A, 3 T 
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Table XXVII. — ;Values of the variable part of the magnetic potential on the 
cosine of the colatitude) and 24 equidistant meridian circles, reckoned from 
is that of the mean summer months. The time is Greenwich noon. 
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Table ^XXVIII. — Values of the variable part of the magnetic potential on the 
cosine of the colatitude) and 24 equidistant meridian circles, reckoned from 
is that corresponding to the mean values for the year. The time is Green- 



JUL = COS u. 


7^ = 0. 


15. 


30. 


45. 


60. 


75. 


90. 


105. 


120. 


135. 


150. 




'98 


- 15^1 


- 7-8 


0^1 


+ 7-1 


+ 13-1 


+ 17-6 


+ 20-5 


+ 21-7 


+ 21-8 


+ 20-8 


+ 19-0 




•94 


~ 44-0 


- 29-3 


12-4 


+ 4^7 


+ 20-1 


+ 32-1 


+ 40-2 


+ 44-4 


+ 45-2 


+ 43-7 


+ 40-7 




•87 


- 99-1 


- 75-5 


- 44-3 


9-8 


+ 23-1 


+ 50-7 


+ 70-1 


+ 80-6 


+ 83-0 


+ 79-8 


+ 72-9 




'77 


-159-6 


128^5 


- 82-7 


-29-3 


+ 23-3 


+ 68-1 


+ 99-8 


+ 116-4 


+ 119-1 


+ 111-7 


+ 99-3 




-64 


-188-2 


-152^6 


99-8 


37-0 


+ 24-0 


+ 74-8 


+ 109-5 


+ 125-9 


+ 126-3 


+ 115-9 


+ 101-3 




•50 


166-2 


130^8 


- 80-9 


-25-8 


+ 25-6 


+ 65-9 


+ 91-0 


+ 100-4 


+ 97-2 


+ 87-8 


+ 77-7 




•34 


110-2 


- 80-9 


44-3 


- 7-2 


+ 23-9 


+ 45-2 


+ 55-5 


+ 56-2 


+ 51-4 


+ 45-7 


+ 43-2 




•17 


- 49-9 


33-3 


14-6 


+ 2-2 


+ 14-7 


+ 21-4 


+ 22-6 


+ 20^2 


+ 16*5 


+ 14-7 


+ 16-1 




•00 


00^0 


00^0 


00-0 


00-0 


00-0 


€0-0 


00-0 


00-0 


00-0 


00-0 


00-0 





VARIATION OF TERRESTRIAL MAGNETISM. 



507 



Earth's surface for 17' latitude circles corresponding to different values of /x (the 
Greenwich towards the East. The period of the year to which the Table refers 
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+ 20-5 


+ 29-9 


+ 38-1 


+ 44^-4 


+ 487 


+ 50-7 


+ 49-7 


+ 45*6 


+ 37-7 




- 15-0 


- 6'8 


+ 1-4 


f 9-6 


+ 17-1 


+ 23-4 


+ 28-3 


+ 31-3 


+ 32-5 


+. 31-6 


+ 28-5 


+ 23*4 


+ 16-6 



Earth^s surface for 17 latitude circles corresponding to different values of [i (the 
Greenwich towards the East. The period of the year to w^hich the Table refers 
wich noon. 





165. 


180. 


195. 


210. 


225. 


240. 


255. 


270. 


285. 


300. 


315. 


330. 


345. 




+ 16*7 


+ 13-5 ' 


+ 9-8 


+ 5-0 


- 0-5 


&'^ 


-13-0 


-18*9 


23-7 


- 26*7 


- 27-4 


- 25*5 


- 21-3 




+ 36*9 


+ 32-2 


+ 26*5 


+ 19*5 


+ 10*1 


1-2 


-14-5 


28*4 


41-6 


- 52-3 


- 58*6 


- 59-6 


- 54-5 




+ 65-2 


+ 57-6 


+ 55-2 


+ 42-2 


+ 31-1 


+ 16*0 


- 4-4 


-28-6 


55-3 


- 80*8 


-100*9 


-112-0 


-111*7 




+ 86-6 


+ 76*4 


+ 69*4 


+ 63*8 


+ 55*0 


+40-4 


+ 16-7 


16*6 


57*2 


- 997 


-137-4 


-163*1 


-171-4 




+ 88-4 


+ 80*6 


+ 78*2 


+ 777 


+ 74-2 


+ 61-8 


+ 36*8 


- 1-9 


51-4 


-104*7 


-1531 


-187*1 


-200-1 




+ 71-8 


+ 71*8 


+ 76*3 


+ 81-3 


+ 80*4 


+ 68-7 


+ 43-4 


+ 3-4 


-45-8 


-97*4 


142-4 


-172*2 


'180-7 




+46*0 


+ 53-6 


+ 63-2 


+ 70*3 


+ 69-9 


+ 58-6 


+ 35*0 


+ 11 


-38-5 


- 77*4 


-108-4 


-125-8 


-126-3 




+ 21*2 


+ 28*7 


+ 36*5 


+41*6 


+41-0 


+ 33*2 


+ 18-6 


- 1-4 


-23*4 


- 43*5 


57-9 


64*1 


- 61-2 




00*0 


00*0 


00*0 


00*0 


00*0 


00-0 


00*0 


00-0 


00*0 


00*0 


00*0 


00-0 


00-0 
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Vii. (Jonciuding Itemarks. 

Faeaday, in the year 1850, discussed tbe diurnal variation of the magnetic needle. 
He showed that the changes which took place during daytime could be accounted for 
by supposing two magnetic poles— namely, a North pole in the Southern hemisphere, 
and a South pole in the Northern hemisphere— to be carried round with the Sun in 
our atmosphere. A glance at fig. 13 will show that our result entirely agrees with 
Faraday^s. The proof that the principal part of the Earth's magnetism is due to 
causes outside its surface would have been almost as complete in the year 1850 as 
it is now, if Faraday had added the remark that, if all three components of the 
variation can be completely accounted for by hypothetical changes taking place outside 
tlie Earth's surface, they cannot be accounted for by changes taking place in the 
inside. 

I cannot agree, however, wdth Faraday in the explanation which he gives of the 
variation. He imagines that the solar radiation, heating up the air, produces a 
sufficient change in its magnetic permeability to account for the observed deflection 
of the lines of magnetic force. 

The magnetic susceptibility of oxygen at the atmospheric pressure and temperature 
is about *5, 10"^, and for air it is smaller still. This would give the magnetic per- 
meability as i '0000006. If the air was entirely removed the change of magnetic 
force would be so small that we could not detect it. I have tried in various wavs to 
find how a partial removal of the atmosphere as a magnetic medium could afiect the 
needle in any appreciable way, but have failed to do so. Faraday suggests that the 
oxygen in the higher regions of the atmosphere might, owing to the greater cold, be 
much more magnetic than what we observe it to be. But, on the other hand, owing 
the smaller density, the permeability would be diminished ; so that I do not think we 
are at present justified in ascribing any material part of the daily variation to a 
change of the magnetic permeability of air due to the heating effect of the Sun. The 
effect of the Moon suggests a tidal action as the cause, arid we may inquire whether 
such a tidal action could produce the observed effects. The late Professor Balfour 
Stewart has suggested that the Earth's magnetic force might induce electric currents 
in the convection currents which flow in the upper regions of the atmosphere. One 
difficulty of this hypothesis was removed by an experimental investigation, by means 
of which I have proved that the air can be thrown into a sensitive state in which 
small electromotive forces will produce sensible electric currents. To bring the air 
into that sensitive state it is only necessary to send an electric current through it 
from some independent source of high potential. It is very likely that the air in the 
upper regions of our atmosphere is in such a sensitive state, and it is quite possible, 
therefore, that the induced electric currents suggested by Professor Balfour Stewart 
really exist. 

The symmetry of the diurnal variation in both hemispheres shows that, if it is due 



510 PEOPESSOR A. SCHUSTER ON THE DIURIS^AL 

to the assumed cause, tbe vertical component of the magnetic force is the important 
factor, as that component changes sign on crossing the magnetic equator. In order 
that electric currents should be induced which could account for the observed move- 
ment of the magnetic needle, it is only necessary to imagine convection currents in the 
upper regions from East to West during certain parts of the day, and from West to 
East at other times. Judging from the analogy of the theory of waves in shallow 
water, a horizontal motion of considerable velocity might be produced by a tidal action 
due to solar and lunar attraction. It is true that no periodic change of the barometer 
has been traced with certainty to a tidal action ; but I suppose that a tidal wave 
must nevertheless exist, and that its horizontal flow might be considerable^ while the 
changes of pressure might escape our attention. As regards the effect of the Sun we 
have, indeed, a daily period of the barometer which is probably due to thermal 
effects. It is curious and suggestive that the horizontal motion which must accom- 
pany the change in pressure is just such as would account for the daily variation of the 
magnetic needle. In the tropics the principal minimum of the barometer takes place 
about 3.40 o'clock in the afternoon, and the principal maximum about 9 o'clock 
in the morning. According to the theory of waves, there would be a horizontal 
movement from West to East in the afternoon, and from East to West in the morning. 
The direction of the induced electric currents would be away from the equator in both 
hemispheres in the afternoon, and towards the equator in the morning. This is 
exactly the system of currents we have been led to, starting from the observed 
magnetic variation. The only difficulty I feel in suggesting that the cause of the 
diurnal variation of the magnetic needle is the diurnal variation of the barometer lies 
in the fact that it would oblige us to place the electric currents into the lower regions 
of the atmosphere, as these only will be much affected by the thermal radiation of the 
Sun. The phase of the barometric oscillation has been found to be reversed on the 
top of mountains, and it would be interesting to see whether the magnetic variations 
show any peculiarities at great heights. "^'^ 

The region of the atmosphere which other considerations lead me to consider as the 
most sensitive to electromotive forces is that of the cirrus clouds, and I should be 
inclined to look to that region for a solution of the question. The lunar action 
seems, according to the researches of Mr. Chambers, to be a modification of the solar 
action rather than an independent effect. This might be accounted for if we suppose 
that the conductivity of the air depends on the position of the Sun, while the electro- 
motive forces depend on the combined positions of the Sun and Moon. 

* [Note added October 11, 1889. — Since writing tlie above I have become acquainted with Hann's 
recent work on tbe diurnal oscillation of the barometer Q Wien, Denkscbriften,' vol. 55, 1889.) It 
appears from the regularity of tbe semidiurnal period in different altitudes and latitudes that its cause must 
lie in atmospberic movements in bigber regions of tbe atmospbere. The reversal of pbase mentioned in 
tbe text is due to local effects and bas notbing to do witb tbe regular oscillation. It seems to be 
exceedingly probable in tbe ligbt of tbese researcbes tbat tbe daily variation of tbe magnetic needle is 
connected witb tbe daily oscillation of tbe barometer in tbe way described in tbe text.] 
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It will be interesting to follow out in future researches the field which this investi- 
gation has opened, especially in order to trace the effect of the Sunspot variation ; but 
for this purpose it is absolutely necessary that different observatories should follow a 
more uniform plan in reducing their observations. It has been found by experience 
that if the hourly readings of the magnetic needle are collected together, and their 
mean taken, that mean is different according as the disturbed days are taken into 
account or rejected. In other words, the disturbances are not irregularly distributed, 
but have a daily period which is mixed up with the regular daily variation. If we 
want to separate the investigation concerning the regular variation and the disturbance 
variation, we must adopt some plan of obtaining the one without the other. I need 
not here describe Sabine's well-known method of doing this. Grave objections have 
been urged against it, but it is still adopted in many observatories. A discussion of 
the various methods of reduction which have been proposed will be found in recent 
Keports of the British Association, and amongst them that adopted by Mr. Wild at 
St. Petersburg seems to me to be the only one which can be justified on strict scientific 
principles. It consists in selecting the curves for the quiet days, of which there are 
always a sufficient number in each month, and not to take account, as Sabine's 
method does, of any reading at all during the disturbed days. We get in this way 
something perfectly definite, namely, the mean variation of the magnetic needle 
during certain specified days. It seems to me that if the heads of different observa- 
tories could adopt some system of intercommunication, by which they could select 
those days which are most quiet all over the world, and if the elements are reduced 
for those days solely at the different stations, we should obtain a series of values for 
different points of the Earth which are strictly comparable with each other. The 
labour of reduction, as far as I can judge, would thereby be seriously diminished. 
The method hitherto adopted at Greenwich is very similar to that of Wild, and will 
not, probably, lead to results which are sensibly different. 

The reduction of the observations by spherical harmonic analysis would be a very 
simple matter according to the method which I have followed, if the results of different 
stations were published in a manner which would lend itself easily to the w^^ork. The 
method of publication adopted at Greenwich is very convenient, and might serve as a 
model to other observatories. Much labour is, however, involved in reducing varia- 
tions in horizontal force and declination to variations in force towards the geographical 
West and North respectively. If all observatories could publish the coefficients of the 
harmonic series of the elements as at Greenwich, but reduced to the geographical 
instead of the magnetic co-ordinates, the progress of magnetic science would be much 
assisted, as every scientific investigation must take the geographical components for 
its starting point. 

I have tried to form an idea as to the degree of accuracy reached in the deter- 
mination of such quantities as the daily variation of declination ; the result is not 
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altogether satisfactory. Mr. Whipple, in the ' British Association Report ' (Birming- 
ham, 1886, page 71), says :— 

^' Contrasting the Kew results with those of Greenwich, we may fairly consider the 
difference to be due in some measure to instrumental causes, the construction of the 
magnetographs being dissimilar at the two observatories. The slight difference in 
position of the two observatories may likewise have some influence.'' 

The difference amounts to about 15 per cent., and it seems as if the question 
whether such a difference can be due to instrumental causes deserves a careful 
examination. Mr. Chambers, at Bombay, has found similarly that the results of the 
magnetographs differ from those obtained by the old magnetometers ; and he seems to 
ascribe the difference to an ^' influence of height above or below the ground level. ^^ 
The height of the magnetometer was 6 feet above ground, and that of the magneto- 
graph 7 1 feet below ; the former gives ranges greater by 7 per cent, for the declination 
variation, and the difference is greater still for the horizontal force component. That 
there should be a real difference of that magnitude in the two positions seems exces- 
sively unlikely, and we must conclude that at present the results given by magneto- 
graphs are doubtful to the extent of about 10 per cent. 

In conclusion, we may sum up the principal results obtained in this paper as 
follows :— 

1. The principal part of the diurnal variation is due to causes outside the Earth's 
surface, and probably to electric currents in our atmosphere. 

2. Currents are induced in the Earth by the diurnal variation which produce a 
sensible effect chiefly in reducing the amplitude of the vertical component and 
increasing the amplitude of the horizontal components. 

3. As regards the currents induced by the diurnal variation, the Earth does not 
behave as a uniformly-conducting sphere, but the upper layers must conduct less 
than the inner layers. 

4. The horizontal movements in the atmosphere which must accompany a tidal 
action of the Sun or Moon or any periodic variation of the barometer such as is 
actually observed would produce electric currents in the atmosphere having magnetic 
effects similar in character to the observed daily variation. 

5. If the variation is actually produced by the suggested cause, the atmosphere 
must be in that sensitive state in which, according to the author's experiments, there 
is no lower limit to the electromotive force producing a current. 

In conclusion, the author begs to return his thanks to Mr. William Ellis for help 
given in some of the calculations, and also to his assistant, Mr. Arthuh Stanton, for 
much labour bestowed on making and checking numerical calculations. 
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Appendix. 

On the Currents Induced in a Spherical Conductor by Variation of an External 

Magnetic Potential 

By Horace Lamb, M,A., FM.S. 

The general formulae for the currents induced in a sphere of uniform conductivity by 
any electric or magnetic disturbances outside it have been given in the 'Phil. Trans./ 
1883, pp. 526 et seq. I here reproduce (with some further developments) so much of 
the investigation as is required for the discussion in Part V. of the foregoing paper. 

Suppose that, the origin being taken at the centre of the Earth, we have an external 
disturbing force, whose magnetic potential near the Earth's surface is 

when O^^ is a solid harmonic of positive degree n, and ?' = y^ •— 1. 

The corresponding values of the components of electric momentum are 

-T^ 1 id d\ 

p 1 ( ^ ^ ^\r\ 

n + 1\ dx 



.J • «' 



yj 1 I d d\ 

~ ^7TT r % ~ '^JxP''' 

the time-factor being omitted here and elsewhere for shortness. For these values 
make 

d¥ dG dM. ^ 

dx dy dz ' 

and also give, for the components of magnetic force, the values 



""■ dy 


da 

dz 


^ --|.W^O,, — (.X 

dn,, 


and, similarly, 




dx ' 

^- dy 

dD,n 
^ dz 
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If II, V, IV be the components of electric current, the equations of induced currents 



are 



pu 



F, 



pv= — G, 



plD 



H, 



where p denotes the specific resistance. Eliminating it, v^ w by means of the relations 



V'F 



iwU, V^Gr 



Anv, 



V^H = — 4:7TW, 



we find that at all points within the sphere F, G, H must satisfy the equations 



(y2 _ p) -p ^ 0, (v' - F) G 



where 



d¥ d^ dK 

clx dy d'Z 



0, 
0, 



4™ . 



(V^ - F) H == 0, 



The appropriate solution of these is 



d 

dz 

d 

dx 

d 



d 
dy 



m 



X 



d 
dz 

d 



] (^n, y 



H-x4H.(^^-;?/,^)^.. 



^ 



where r = {x^ + 2/^ + ^^)% ^n denotes a solid harmonic of degree n^ and 



Xn (0 = 1 + ^ 



-.+ 



r 



4 



2?^ + 3 • 2.4.27^ + 3.2^ + 5 

d Y sinh f 



"T" . . « 



3.5 .. . 2/2 4- 1. 



?^r/ 



The total magnetic potential outside the sphere will be 

where O^^n^i is the part due to the indu.ced currents. The values of F, G, H at 
external points will then be 



F 



G 



II 



y 



d 



d 



d: 



z ■ 



'y. 



/• d 
\z 



X 



dx 

d 
dy 



X 



dy, 

dz^ 
d 



12 -' 



y 



dx 



n -{- 1 
1 

1 

n + 1 ^^-^ 



1 

n 



u—i h 
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k) Jl tJ 



It remains to introduce the conditions to be satisfied at the surface of the sphere 

• • • 

(r = R). The continuity of electromotive force, i.e., of ¥, G, H, requires 



Xn {^^) 



0) 



n 



9^ + 1 



n 



0_;^_i. [r = R] 



The continuity of the magnetic force involves the continuity of the space- 
derivatives of F, G, H, and, therefore, of dF/dr, dG/dr, dH/dr, Hence 

{k^x,, (kU) + o%Xn {kB) } 0)^ = ;;7^~ O,, + --^^ fl«;,„i . [r = R] 



n -i- 1 



n 



thence find 



k~R . x^ (kU) . 0)^ = — ;^ — 0_,,_i, 



n 



{kux'n im + {2n + 1) xn {m ] CO,, = -^ a„, 



V 



n -i- 1 



which are equivalent to 






7rE 



•^ 



2n + 1.2?^ + 3 



X^ + l (^-tt) ' ^n — ^ 0_;^_i , 



V 



Xn^l (/^R). 0)^ 



n + 1 



•O.. 



7* 



Hence 



ix _K— 1 



k^R' 



O 



ti 



%'M + l(M) 



r =: R] 



R 



This gives the ratio of surface potentials, and, therefore, of horizontal forces, due 
to internal and external iufluences respectively. Since this ratio is '^complex," there 
will be a difference of phase, as well as of amplitude. The corresponding ratio of 
vertical forces is 



(a/'aL __ M _ 1 



dr 



da 



n 



dr 



^ + 1 O 



n — i 



n 



n 



n 



2n + 1.271 + 3* %«^i(M) 



r = R 



To interpret these results it is necessary to calculate the function 



r 



where 



2^^+1.27^ + 3 Xn-l(0 



2 



¥R 



a 



47r|?E^ . 



p 

3 IT 2 



i = iB, say. 
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For moderate values of 8 we may use the form 



where 






n 



^8 



1 - 



W 



A 4- I'R 

271 + 1.2w + 3 A„_i + i'^n-\ 



+ 



•o 

XJ^ - — 



2.4.2^1 + 3.2^ + 5 ' 2.4.6.8/2n, + 3.2^2, + 5.2^ + 7.2fi + 9 



2.27^ + 3 2.4.6.2% + 3.29^ + 5.27^ + 7 



* • « * 



The following Table gives the values of A^, B^, A3, B3, A5, Bg, for various values of 8. 
It may possibly be of service in other investigations. 



lABLE JLA.iJv. 



^, 


Ai. 


Bp 


A3. 


Ba. 


A5. 


B5. 


1 


+ 


0-996429 


"- 


0-099934 


+ 0-998738 


+ 0-056639 


+ 0-999369 


+ -038466 


2 




0-985726 




0-199471 


+ 0-994952 


+ 0-110982 






3 




0-967918 


— 


0-298216 


+ 0-988647 


+ 0-116230 






4 


+ 


0-943060 




0-395773 


+ 0-979832 


+ 0-221187 




\ 


5 


+ 


0-911184 




0-491751 


+ 0-968618 


+ 0-275767 






6 


+ 


0-872401 


+ 


0-585759 


+ 0'954720 


+ 0-329843 






7 


4- 


0-826800 




0-677412 


+ 0-938455 


+ 0-383350 






8 




0-774498 




0-766327 


+ 0-919714 


+ 0-436182 






9 


+ 


0-715628 




0-852126 


-f 0-898610 


+ 0-488246 






10 




0-660341 


_-. 


0-934439 


+ 0-875083 


+ 0-539448 


+ 3-936309 


+ -3784 


20 




0-310366 




1-489206 


+ 0-516311 


+ 0-98^a36 






30 




1-628644 


+ 


1-351846 


-- 0-029614 


+ 1-248628 






40 


- 


2-912855 




0-337304 


0-688914 


+ 1-266555 






60 




3-716064 





1-663371 


- 1-366602 


+ 0-993199 






100 


+ 


12-514840 




10-691120 


-^ 1-308236 


3-811231 


- 1-9434 


-m^^^ 



If we write 



the above fraction becomes 



^n \ ^ ^n — ^rfi "j 






2?^ + l,27^ + 3 S,,„i 

If we prefix the minus sign this gives the ratio (^') of the vertical forces, 
the ratio (^) of the horizontal forces we must multiply by ni{n + 1). 
For large values of 8 we may make use of the second form of x#- Thus 



To get 



Xi 



Xs 



cosh ^ sinh ^ 



^3 

15\ 




>■> 



^4 + ^ ! cosh t 



6 



. ^. + "i" ) sinh t 



whence 
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^' %3 (D ^ (1 + wt') cosh r- (6rv+ isr^) sinh r 

5 • 7 %i (?) cosh ? — ?-i sinh f 



Here 



^ = {iSf = (1 + i) SVx/2 = (1 + i) /3, say. 



If yS is moderately large we may put cosh 1/ sinh ^ = 1 approximately. The error 
thus committed is of the order e""^^ Since the value of e""-^^ has six cyphers after the 
decimal point, this approximation is amply sufficient for /3 > 7, or say for § > 100. 
The above fraction is then 

^ i-r~^ ^ (1-1^-1) + ^4/3-^ 

It is by these methods that Tables XX. and XXI. above were calculated. The 
values of p, the specific conductivity, given in the fifth columns, were obtained from 
the formula 

by putting 

27rR = 410^ cm., 2it/p = 86,400/m sees., 

where m denotes the number of complete periods in a day, and is therefore = 1 for 
the diurnal and = 2 for the semidiurnal variations. 

As the resistance diminishes, the diflference of phase tends to zero, and the ratio 
of normal forces to the value •— 1 ; i,e,, the total normal force at the surface tends to 
zero, in accordance with the theory of electromagnetic screens. 

The ratio of the total vertical to the total horizontal force in any assigned 
direction is 

d 
ar 



d 



where ^drj denotes a linear element drawn in the proper direction on the Earth's 
surface. By means of the preceding results this can be put in the form 

n (71 4- 1) Xn (H) . / diln 

The coefficient may be calculated independently, by a proper adaptation of the 
previous methods, or we may deduce its value from the results already obtained, in 



518 ON THE DIURNAL VARIATION OF TERRESTRIAL MAGNETISM. 

the mannei^ explained by Professor Schuster. The fanction actually computed by 
Professor Schuster in Tables XXII. and XXIII. is the ratio 



tv^Li^, 



dr 
which 



/^R;^/(/^R) + {n^l)xn(]c^) 



For large values of 8, i.e. for sufficiently small values of p, we may put 
cosh ^ = sinh ^, whence, keeping only the most important term in ^n (0? ^^^ fraction 
written becomes 

= (7^ + l)/C = ^•e--/*. 

The difference of phase therefore tends to the limit 45°, as remarked by Professor 
Schuster. For 8 = 100, this formula gives for the reduction of amplitude the value 
•3 and '5 in the cases n = 2 and n == 4 respectively. 



